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HIGHLIGHT  SUMMARY 


Drilled  piers  or  bored  piles  can  be  an  efficient 
means  to  control  slope  movements  where  other  corrective 
measures  fail  to  insure  stability  or  when  their  use  is 
prohibited  due  to  space  limitations.   A  design  metho- 
dology is  proposed  herein  to  assess  both  the  forces 
acting  on  the  piles  (or  piers)  and  the  influence  of  the 
pile  row  on  the  stability  of  the  slope.   The  forces 
acting  on  the  pile  sections  above  and  below  the  criti- 
cal surface  are  calculated  using  the  theories  of  plas- 
tic deformation  and  subgrade  reaction,  respectively. 

The  theory  of  plastic  deformation  is  based  on  the 
assumption  that  the  soil  surrounding  the  piles  is  in  a 
state  of  plastic  equilibrium.   Under  this  assumption, 
the  force  acting  on  the  passive  piles  is  expressed  as  a 
function  of  the  soil  strength  parameters  and  of  the 
pile  diameter,  spacing  and  position.   The  actual  reac- 
tion force  exerted  by  the  piles  on  the  slope  is  assumed 
to  be  a  fraction  of  the  force  corresponding  to  the 
plastic  state  condition.   The  Friction  Circle  Method  of 


xvi 


slope  stability  analysis  is  extended  to  incorporate  the 
reaction  force  and  determine  the  critical  surface  and 
safety  factor  of  the  reinforced  slope.  A  computer  pro- 
gram is  also  developed  to  obtain  the  displacement, 
bending  moment,  and  shear  force  profiles  along  the 
piles.  A  step-by-step  procedure  is  proposed  to  achieve 
an  optimum  design  solution  which  provides  a  required 
factor  of  safety  for  the  slope  and  insures  the 
integrity  of  the  piles. 


CHAPTER  I.   INTRODUCTION 


Insuring  the  stability  of  both  natural  and  man- 
made  slopes  continues  to  be  an  important  problem  in 
geotechnical  engineering.   Slides  account  for  many 
civil  engineering  failures  and  often  result  in  exten- 
sive property  damage  and  sometimes  loss  of  human  life. 
There  is  no  universally  accepted  method  for  the  preven- 
tion and/or  correction  of  landslides.  Each  slide  is 
unique  and  should  be  considered  on  the  basis  of  its  own 
individual  characteristics.  Avoidance  of  a  potential 
slide  area  can  be  a  primary  consideration  when  select- 
ing a  new  site.   Otherwise,  corrective  measures  can  be 
taken  which  include:   (1)  improving  the  slope  geometry 
by  changing  the  slope  angle,  excavating  the  soil  at  the 
head,  or  increasing  the  load  at  the  toe;  (2)  construct- 
ing a  compound  slope;  and  (3)  providing  surface  and 
subsurface  drainage.  However,  where  such  corrective 
measures  fail  to  insure  stability  or  when  their  use  is 
prohibited  due  to  space  limitations,  retaining  struc- 
tures may  be  necessary.  Piles  carried  across  the 
active  or  potential  failure  surface  can  be  used 


efficiently  as  the  retaining  element,  since  they  can 
often  be  installed  without  disturbing  significantly  the 
equilibrium  of  the  slope. 

The  analysis  of  a  slope/pile  system  requires  that 
the  soil  pressure  acting  on  the  piles,  and  the  subse- 
quent reaction  that  the  piles  provide  to  the  slope,  be 
known.   There  are  several  techniques  which  are 
currently  used  to  estimate  the  pressures  that  act  on 
active  or  passive  piles.   They  include  the  coefficient 
of  subgrade  reaction  method,  methods  based  on  the 
assumption  of  elastic  soils  and  other  empirical 
methods.  Most  of  these  techniques  either  consider  the 
solution  of  a  single  pile  and  introduce  corrective  fac- 
tors to  analyze  the  row  of  piles,  or  treat  the  row  of 
piles  as  a  retaining  wall.   The  first  approach  can  lead 
to  inaccurate  answers  in  case  of  passive  piles,  while 
the  second  approach  provides  very  conservative  answers. 
Moreover,  very  limited  research  has  been  done  to  deter- 
mine the  growth  mechanism  of  the  lateral  force  acting 
on  a  pile  embedded  in  a  slope  and  the  increase  in  the 
factor  of  safety  of  the  slope  due  to  the  presence  of 
the  piles.   In  addition,  a  design  methodology  is 
required  to  assist  in  selecting  design  parameters,  such 
as  pile  diameter,  spacing,  and  location  that  will 
insure  the  stability  of  both  the  slope  and  the  pile. 


In  this  report,  the  pressure  acting  on  the  piles 
is  calculated  using  the  theory  of  plastic  deformation 
developed  by  Ito  and  Matsui  (1975).  The  major  assump- 
tion in  this  theory  is  that  the  soil  in  the  area  sur- 
rounding the  piles  is  in  a  state  of  plastic  equili- 
brium. With  this  assumption,  the  force  acting  on  the 
piles  can  be  expressed  as  a  function  of  the  soil 
strength,  and  the  pile  diameter,  spacing  and  location. 
The  safety  factor  of  the  slope  after  the  placement  of 
the  piles  is  calculated  assuming  that  a  portion  of  that 
force  is  counteracting  the  driving  forces  of  the  slope. 
Finally,  a  step-by-step  procedure  is  proposed  to  select 
design  parameters,  which  will  achieve  the  stability  of 
the  slope  and  adequate  dimensioning  of  the  piles. 

The  concepts  and  methods  presented  herein  are 
incorporated  in  two  computer  programs,  which  should 
provide  efficient  design  tools  to  practicing  engineers. 
The  first  computer  program  calculates  the  safety  factor 
of  the  reinforced  slope  as  a  function  of  pile  diameter, 
location  and  spacing.  The  second  program  calculates 
the  shear  force,  bending  moment  and  displacement  along 
the  length  of  the  pile  as  a  function  of  the  soil  pres- 
sure, and  the  pile  diameter,  stiffness,  location,  spac- 
ing, and  boundary  conditions.  These  programs  can  be 
used  iteratlvely  to  provide  an  optimum  design  solution. 


CHAPTER  II.   SLOPE  REINFORCEMENT 


Soil  reinforcement  techniques  such  as  reinforced 
earth,  stone  columns,  soil  anchors,  and  piles  and  piers 
can  be  used  to  stabilize  slopes  and  highway  embank- 
ments.  When  properly  designed  and  constructed,  rein- 
forcing elements  can  be  appropriate  to  solve  stability 
problems,  particularly  where  space  Is  restricted,  and 
they  can  provide  cost-effective  solutions  to  many  tran- 
sportation design  and  construction  problems. 

Piles,  carried  across  an  active  or  potential 
failure  surface  and  aligned  to  form  a  cantilever  wall, 
have  been  used  in  the  past,  although  early  uses  of  such 
walls  have  resulted  in  some  failures  (Baker  and 
Marshall, 1958;  Root, 1958;  Broms,1969).  Merriam  (1960) 
described  the  unsuccessful  attempt  to  stabilize  the 
Portuguese  Bend  slide  of  the  Palos  Verdes  peninsula  in 
Los  Angeles  County,  where  an  ancient  slide  mass  had 
been  reactivated  by  a  new  housing  development  in  the 
area.  Movements  at  the  rate  of  1.25  in.  per  day  were 
observed.   In  order  to  stabilize  the  slide,  25  4-ft 
diameter  concrete  caissons  were  installed  in  separate 


groups  of  2  or  3.   Initially  the  soil  flowed  around  the 
caissons,  but  these  were  eventually  broken  by  the  mass 
movement.  Since  the  original  slide  had  remolded  the 
soil  to  its  residual  strength,  the  small  relative 
increase  in  shear  strength  provided  by  the  caissons  had 
an  insignificant  effect  on  the  slope  stability. 

Since  then,  more  effective  walls  have  been  con- 
structed.  Shallow  slides,  up  to  20-ft  deep,  have  been 
stabilized  by  the  use  of  sheet  piling  (Toms  and 
Bartlett,  1962)  and  closely  spaced  driven  cantilever 
piles  (Zaruba  and  Mencl,1967).  A  root  pile  wall  was 
used  to  correct  a  landslide  near  Monessen,  Pennsylvania 
(Dash  and  Jovino, 1980).   Stone  columns  have  been  used 
effectively  in  Europe  (Goughnour  and  DiMaggio, 1978). 
Drilled  cantilever  pier  walls  have  been  used  success- 
fully in  the  Ohio  River  Valley  (Nethero,1982).   Satis- 
factory results  were  achieved  with  18-in.  to  30-in. 
diameter  piers  spaced  5  to  7-ft  on  center.  Actual  instal- 
lations of  such  piers  include  the  use  of  cast-in-place 
concrete  cylindrical  columns  in  an  effort  to  stop  soil 
movements,  which  had  been  sufficient  to  cause  closing 
of  a  roadway.   In  a  different  application,  piers  placed 
at  the  toe  of  a  slope  stopped  a  slide  which  had  been 
triggered  by  an  unbraced  cut. 

Deep-seated  slides  have  also  been  stabilized  suc- 
cessfully by  large  diameter  cylinder  piles,  anchored 


sheet  pile  walls,  or  rock  anchored  cylinder  piles.  Some 
of  these  cases  will  be  discussed  below. 

Large  diameter  cast-in-place  concrete  cylinders 
were  used  for  a  sidehill  stabilization  during  the  con- 
struction of  the  Seattle  freeway  in  Washington  (Andrews 
and  Klassel, 1964).   During  excavation  for  the  footing 
of  a  bridge  structure,  a  large  deep-seated  slide 
started  to  develop,  immediately  adjacent  to  a  seven- 
story  apartment  building.  A  retaining  wall,  which  had 
been  planned  between  the  bridge  footing  excavation  and 
the  building,  could  not  be  built  because  further  exca- 
vation for  the  wall  would  cause  the  slide  to  progress 
and  endanger  the  building.   To  avoid  this,  cast-in- 
place  cylinder  piles  were  used  to  form  the  wall.  To 
develop  the  resistance  required  to  withstand  the 
assumed  loads,  the  design  called  for  4.75  ft  cylinders 
on  6-ft  centers.   The  cylinders  were  designed  to 
penetrate  the  slip  plane  and  other  potential  sliding 
planes.  Welded  steel  beams  were  employed  to  reinforce 
the  concrete  cylinders. 

A  sheet  pile  wall  anchored  with  pre-tensioned  soil 
anchors  helped  to  stabilize  a  potentially  dangerous 
slope  during  construction  in  the  Ohio  River  valley. 
According  to  D'Appolonia  et  al.  (1967),  the  wall  pro- 
vided short  term  stability  by  preventing  progressive 
failure,  while  drainage  assured  long-term  stability. 
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Cylinder  piles  were  used  in  an  effort  to  stop 
movement  during  construction  for  the  freeway  on  Potrero 
Hill  in  San  Francisco  (Nicoletti  and  Keith, 1969).   A 
rock  bolted  retaining  wall  was  installed  to  restrain  an 
unbraced  cut.   The  top  of  a  railroad  tunnel  was  located 
25-ft  below  the  freeway  and  near  the  top  of  the  wall. 
The  overall  stability  of  the  tunnel,  wall,  and  the 
slope  rising  above  it,  had  been  unbalanced  by  the  cut. 
Towards  the  end  of  the  construction,  large  lateral 
movements  were  observed  in  the  tunnel.   In  order  to 
stop  the  movements,  a  wall  was  constructed  of  a  series 
of  heavy  steel  piles,  placed  in  drilled  holes  on  each 
side  of  the  tunnel,  and  connected  by  steel  struts 
across  the  top  of  the  tunnel.  This  technique  proved  to 
be  efficient  in  preventing  further  movement. 

A  double  row  of  anchored  piles  became  necessary 
during  railroad  work  in  Belgium  (DeBeer  and  Wal- 
lays,1970).   In  order  to  enlarge  the  railroad  bed,  an 
existing  slope,  nearly  at  its  equilibrium,  was  cut 
back,  maintaining  the  original  angle.  The  excavations 
reactivated  slips  which  had  previously  occurred  in  a 
crushed  schist  mass,  and  a  series  of  retrogressive 
slips  was  initiated.  Since  the  slope  was  located  in  a 
city,  only  limited  flattening  was  possible. 
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Hence,  the  stability  was  insured  by  the  resistance  pro- 
vided by  3.5  to  5-ft  diameter  piles  spaced  6.5  to  10-ft 
on  centers. 

An  anchored  cylinder  pile  wall  was  installed  dur- 
ing construction  of  interchanges  to  connect  1-471  with 
local  traffic  arteries  in  Cincinnati,  Ohio  (Offen- 
berger,  1981).  Excavations  on  the  southwest  side  of 
Mt.  Adams  triggered  large  ground  movements  in  an  area 
known  to  have  experienced  small  movements  in  the  past. 
Data  from  slope  indicators  showed  that  the  earth  move- 
ment was  in  a  deep  seated  weak  clay  layer  near  the  rock 
surface,  and  well  below  the  elevation  of  retaining  wall 
footings.  A  cantilever  wall  was  determined  to  be  inap- 
propriate, because  considerable  embedment  into  the  rock 
would  be  necessary.  In  addition,  undesirable  downslope 
displacement  would  be  needed  to  develop  the  required 
lateral  resistance.   Thus,  a  cylinder  pile  wall  with 
tie-back  tendons  was  recommended.   The  cylinder  piles 
were  drilled  through  50  ft  of  overburden  and  socketed 
into  shale  and  limestone. 

Recent  applications  of  slope  stabilization  using 
steel  piles  have  been  described  by  Ito  and  Matsui 
(1981).  Several  landslides  have  been  stabilized  in 
Japan  using  such  piles.   In  the  United  States,  piers 
with  tie-back  rock  anchors  were  drilled  through  a 
landslide  and  embedded  in  an  unweathered  schist  during 


construction  at  the  Geyser  Power  Plant  in  California 
(Hovland  et  al. ,  1982). 

These  case  histories  indicate  that  the  use  of 
piles  for  stabilizing  slopes  is  becoming  more  common. 
However,      little  information  is  currently  available 
regarding  (1)  the  pile  behavior  under  lateral  loading 
induced  by  the  movement  of  the  slide  and  (2)  the  effect 
of  the  pile  row  on  the  overall  stability  of  the  rein- 
forced slope.   The  research  reported  herein  is  a  step 
towards  solving  these  two  problems  and  developing  a 
methodology  for  the  analysis  and  design  of  the 
slope/pile  system. 
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CHAPTER  III.   SLOPE  STABILIZATION  USING  PILES 


The  analysis  of  a  slope  reinforced  with  piles 
requires  that  the  force  applied  to  the  piles  by  the 
failing  mass,  and  in  turn  the  reaction  from  the  piles 
to  the  slope,  be  known.   In  addition,  a  slope  stability 
analysis  that  takes  into  account  the  reaction  force  is 
necessary.  In  this  chapter,  a  theoretical  method  for 
the  calculation  of  pressures  acting  on  passive  piles 
(piles  subjected  to  lateral  loading  by  horizontal  move- 
ments of  the  surrounding  soil)  is  discussed.  Then,  the 
resulting  pile  reaction  is  incorporated  into  a  slope 
stability  analysis  to  determine  the  factor  of  safety  of 
a  slope  reinforced  with  a  row  of  piles.   Finally, 
parametric  studies  are  performed  to  assess  the  influ- 
ence of  the  location  and  spacing  of  the  piles  on  the 
stability  of  the  slope. 
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PILES  IN  SOIL  UNDERGOING  LATERAL  MOVEMENT 


The  problems  of  active  piles  (piles  subjected  to  a 
horizontal  load  at  the  head  and  transmitting  this  load 
to  the  soil)  embedded  in  clay  or  sand  and  subjected  to 
horizontal  static  loads  of  short  or  long  duration  have 
been  treated  by  several  authors  (Matlock  and  Reese, 
1960;  Davisson  and  Gill,  1963;  Broms,  1964;  Vesic, 
1965;  Davisson,  1970;  Poulos,  1973;  Jamiolkowski  and 
Garassino,  1977;  Kishida  and  Nakai,  1977;  Coyle  et  al. , 
1983).   By  making  simplifying  assumptions  with  regard 
to  the  deformation  characteristics  of  the  soil  layer 
surrounding  the  piles,  these  authors  arrived  at  accept- 
able solutions  of  piles  subjected  to  horizontal  loads 
at  the  top.   In  most  cases,  the  problem  of  a  single 
pile  has  been  solved,  and  some  correction  factors  for 
group  effects  have  been  introduced. 

For  the  case  of  passive  piles  (piles  subjected  to 
lateral  soil  movements),  the  problems  are  more  compli- 
cated because  the  lateral  forces  acting  on  the  piles 
are  now  dependent  on  the  soil  movements,  and  these  are 
in  turn  affected  by  the  presence  of  the  piles.  For 
example,  it  is  possible  for  a  pile  group  to  stop  these 
movements,  creating  a  completely  different  situation 
than  that  for  a  single  pile.  Hence,  the  solution  for  a 
single  pile  cannot  be  easily  adapted  for  the  situation 
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of  a  pile  group,  although  several  authors  have  sug- 
gested such  an  approach  (Poulos,1973;  Baguelln  et 
al.  ,1976;  Viggiani.,1981).  Other  researchers  have  con- 
sidered the  problem  from  the  fundamental  standpoint  of 
group  (row)  action  (Ito  and  Matsui,  1975;  Winter  et 
al.,  1983).  Several  theoretical  approaches  for  the  cal- 
culation of  pressures  on  piles  placed  in  deforming 
soils  will  be  presented  here. 

Theoretical  Solutions 

Analytical  models  that  have  been  used  to  obtain  a 
theoretical  solution  for  piles  placed  in  deforming 
soils  include:  1)  a  soil  characterized  by  a  modulus  of 
subgrade  reaction;  2)  a  soil  considered  to  be  an 
elastic-plastic  material;  3)  the  finite  element  method 
used  with  bilinear  or  hyperbolic  approximations  of  the 
soil  behavior;  and  4)  other  approaches. 

Subgrade  Reaction  Method 

In  the  modulus  of  subgrade  reaction  method,  the 
pile  is  considered  to  be  a  beam  on  an  elastic  founda- 
tion, and  the  equilibrium  equation  is  written  as: 


dz         r 
where    EI  ■  pile  stiffness 
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v  =°  horizontal  displacement  of  the  pile  at 

P 

depth  z 

v  =»  horizontal  displacement  the  soil 
's 

would  undergo  without  the  pile 

K  ■  subgrade  reaction  modulus,  function  of  z 
s 

and  (yp  -  yg) 


Computer  programs  are  available  to  solve  this 
equation  (Baguelin  et.al.,  1976),  but  problems  exist  in 
the  evaluation  of  y  .  When  measured  values  of  yg  are 
introduced,  the  solution  of  this  equation  agrees  well 
with  experimental  results.  When  theoretical  values  of 

v  are  introduced,  the  solution  of  the  equation  is 

1  s 

unrealistic  (Theoretical  values  of  yg  are  discussed  by 
Poulos,  1967  and  Canizo  and  Merino,  1977;   Tables  can 
be  used  to  determine  the  horizontal  displacements  of  an 
elastic  layer  subjected  to  various  vertical  load  condi- 
tions). In  addition,  the  applicability  of  the  solution 
to  a  row  of  piles  is  uncertain  because  the  group  action 
modifies  the  initial  soil  deformation  values  that  are 
used  in  the  calculations  (DeBeer,  1977). 


Other  researchers  have  used  the  subgrade  reaction 
concept  for  the  calculation  of  the  horizontal  force. 
Fukuoka  (1977)  studied  the  pile-soil  interaction  in 
slopes  subject  to  creep  by  means  of  a  Winkler  type 
model.  He  assumed  that  the  unbalanced  force  that 
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causes  the  creep  is  a  function  of  the  velocity  of  the 
creeping  mass.  Measurements  are  needed  along  the 
length  of  the  pile  for  the  determination  of  such  a 
velocity,  and  in  turn,  for  the  calculation  of  the 
unbalanced  force.  A  reaction  force  provided  by  the 
piles  and  equal  to  the  unbalanced  force  should  be  ade- 
quate to  stop  the  landslide.  To  calculate  the  lateral 
resistance  of  the  piles  against  the  landslide,  Fukuoka 
utilized  the  modulus  of  subgrade  reaction  method,  and 
gave  equations  for  the  calculation  of  the  deflection 
curve  of  the  pile  and  for  the  bending  moments  and  hor- 
izontal reactions.  However,  when  the  soil  is  moving, 
it  is  very  difficult  to  define  a  coefficient  of 
subgrade  reaction  to  be  used  in  the  equations  (DeBeer, 
1977). 

Viggiani  (1981)  used  the  subgrade  reaction  theory 
to  analyze  the  interaction  between  a  sliding  mass  and  a 
stabilizing  pile.   His  approach  is  based  on  concepts 
developed  by  Broms  (1964)  who  evaluated  the  ultimate 
capacity  of  a  vertical  pile  acted  upon  by  a  horizontal 
load,  using  an  estimated  subgrade  reaction  coefficient. 
Viggiani  evaluated  the  maximum  shear  force  that  a  pile 
can  receive  from  the  sliding  mass  and  transmit  to  the 
underlying  soil,  by  analyzing  six  possible  failure 
mechanisms  for  the  pile.   This  theory  can  be  used  for 
the  design  of  the  piles  against  an  ultimate  pressure. 
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However,  it  does  not  provide  a  way  to  calculate  a  value 
of  the  pressure  at  the  initial  stage  of  the  landslide, 
which  could  be  used  to  determine  the  effect  of  the 
piles  on  the  factor  of  safety  of  the  slope. 

Elastic-Plastic  Material 

Calculations  based  on  the  assumption  of  an 
elastic-plastic  material  are  described  by  Poulos 
(1973).   The  solution  of  the  problem  is  obtained  by 
imposing  compatibility  of  displacements  between  the 
pile  and  the  adjacent  soil.   The  pile  displacement  is 
evaluated  from  the  bending  equation  of  a  thin  strip. 
The  soil  displacement  is  evaluated  from  the  Mindlin 
equation  (Mindlin,  1936)  for  horizontal  displacement 
caused  by  horizontal  loads  within  a  semi-infinite  mass. 
By  solving  a  system  of  equations,  the  displacement  and 
horizontal  pressure  on  the  pile  can  be  evaluated.  This 
method  is  applicable  only  when  dealing  with  an  ideal 
soil  and  a  single  pile. 

More  recently,  Oteo  (1977)  refers  to  a  method 
applied  by  Begemann-DeLeeuw  (1972).  The  method  has 
been  derived  for  the  analysis  of  pressures  on  piles 
when  the  surrounding  soil  undergoes  horizontal  dis- 
placement due  to  a  surficial  surcharge  in  the  vicinity. 
It  considers  the  soil  to  be  a  linearly  elastic  material 
and  is  based  on  stress  and  displacement  fields  given  by 
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the  Boussinesq  theory.   This  technique  is  only  applica- 
ble to  the  case  of  a  known  surcharge  applied  to  the 
slope. 

Finite  Element  Method 

The  applicability   of   finite-element   techniques    to 
the  analysis   of   slope  stabilization  has   been  discussed 
by  Rowe  and  Poulos    (1979).      They   employed  a  two- 
dimensional   finite   element   model   for   soil-structure 
interaction  proposed  earlier  by  Rowe  et  al.(1978).      It 
is    recognized   that   a  finite   element   analysis   of    the 
stabilization  problem  should  make  allowance  for  the 
soil-structure   interaction  effects,    but   also   for 
three-dimensional  effects  such  as   arching  between 
piers.     A  computer  program  which  can  model   these 
effects    is   being  developed  at   Purdue  University    (Oak- 
land and  Chameau,    1983).      The  program  can  also  evaluate 
the  influence   of   piles   position,    size,    spacing  and 
stiffness    on  slope   movements.      When   completed,    this 
program  will  be  added  to  the  design  methodology  pro- 
posed herein  to  form  a   complete   design-analysis  pack- 
age. 

Other  Approaches 

Several  other  methods  have  also  been  used  in  the 
past.   Baker  and  Yonder  (1958)  suggested  calculating 
the  thrust  against  the  piles  by  the  procedure  of  slices 
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and  considering  the  piles  as  cantilever  beams,  provided 
that  they  penetrate  into  a  stable  soil  layer  for  one 
third  of  their  total  length.  Andrews  and  Klassel 
(1964)  and  Nethero  (1982)  calculated  the  pressure  on 
the  pile  by  using  the  concept  of  passive  and  active 
pressures  assuming  that  the  piles  form  a  cantilever 
wall.  However,  analyzing  the  pile  group  as  a  retaining 
wall  can  lead  to  very  conservative  designs,  since  the 
soil  arching  between  the  piles  is  not  taken  into 
account. 

DeBeer  and  Wallays  (1970)  reported  a  method  for 
determining  the  forces  and  bending  moments  on  a  pile  in 
the  case  of  an  unsymmetrical  surcharge  placed  around 
the  pile.   The  method  is  based  on  concepts  introduced 
by  Brinch  Hansen  (1961)  on  the  ultimate  resistance  of 
rigid  piles  against  transverse  forces. 

Wang  et.al.  (1979)  proposed  a  semi-empirical  tech- 
nique to  calculate  the  pressure  on  piles  embedded  in 
deforming  soils  using  inclinometers  to  determine  the 
soil  movements.  The  soil  pressures  required  to  cause 
those  movements  where  then  estimated  using  the  modulus 
of  subgrade  reaction  theory. 

Winter  et.al. (1983)  described  a  method  for  slope 
stabilization  using  piles  based  on  the  viscosity  law  of 
cohesive  soils.   It  uses  a  solution  derived  for  the 
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horizontal  pressure  against  piles  in  viscous  soils. 
Two  requirements,  a  desired  reduction  of  the  sliding 
velocity  and  a  safe  maximum  bending  moment  on  the 
piles,  are  satisfied  by  the  correct  choice  of  pile 
spacing  and  pile  diameter.  This  method  is  only  appli- 
cable for  the  case  of  a  clay  mass  subject  to  creep. 

The  Theory  of  Plastic  Deformation 

In  most  of  the  methods  discussed  so  far,  a 
theoretical  solution  for  the  passive  piles  was  obtained 
by  either  solving  the  problem  of  a  single  pile  or  by 
analyzing  the  pile  group  as  a  retaining  wall.  The 
method  proposed  by  Winter  et.al.(1983)  considers  the 
solution  of  piles  placed  in  a  row  and  takes  into 
account  the  spacing  between  the  piles  at  the  beginning 
of  the  analysis.  However,  this  method  can  only  be  used 
in  purely  cohesive  slopes  undergoing  creep.  A  theoreti- 
cal method  has  been  proposed  by  Ito  and  Matsui  (1975) 
to  analyze  the  growth  mechanism  of  lateral  forces  act- 
ing on  stabilizing  piles  when  the  soil  is  forced  to 
squeeze  between  the  piles.  In  this  analysis,  the  plas- 
tic state  which  satisfies  the  Mohr-Coulomb  yield  cri- 
terion is  assumed  to  occur  in  the  soil  just  around  the 
piles.  The  method  was  developed  specifically  to  calcu- 
late pressures  that  act  on  passive  piles  in  a  row.  The 
assumptions  of  the  theory  of  plastic  deformation  are 


19 


given  in  the  next  section.   Then,  a  derivation  of  its 
application  to  piles  used  in  slope  stabilization  is 
presented.  Finally,  the  validity  of  the  theory  is  dis- 
cussed. 

Assumptions 

It  is  assumed  that  the  piles  placed  in  plastically 
deforming  ground,  such  as  a  landslide,  can  prevent 
further  plastic  deformations.   In  order  to  design  the 
piles,  the  lateral  forces  need  to  be  estimated  as  accu- 
rately as  possible.   These  forces,  however,  are  a  func- 
tion of  the  movement  of  the  sliding  mass.   They  may 
vary  from  zero,  in  case  of  no  movement,  to  an  ultimate 
value,  in  case  of  large  movements.  Ito  and  Matsui 
(1975)  developed  the  theory  of  plastic  deformation  to 
estimate  a  value  of  the  lateral  force  between  these  two 
extremes,  assuming  that  no  reduction  in  shear  resis- 
tance along  the  sliding  surface  has  taken  place  due  to 
strain-softening  caused  by  the  movement  of  the 
landslide.  For  that  reason,  only  the  soil  around  the 
piles  (Figure  1)  is  assumed  to  be  in  a  state  of  plastic 
equilibrium,  satisfying  the  Mohr-Coulomb  yield  cri- 
terion. Hence,  the  lateral  load  acting  on  the  piles 
can  be  estimated  regardless  of  the  state  of  equilibrium 
of  the  slope. 

The  theory  of  plastic  deformation  is  based  on  the 
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F     -r- 


Direction  of 
deformation 


Figure    I      Plastically  Deforming  Ground  Around  Stabilizing 
Piles  (After  Ito  and  Matsui,    1975). 
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following  assumptions: 

1)  When  the  soil   layer  deforms,    two  sliding  sur- 
faces,  AEB  and  A'E'B',    occur,    making  an  angle   of 
(tt/4  +    <J»/2)   with   the   x-axis    (Figure    1); 

2)  The  soil  is   in  a  state   of  plastic  equilibrium 
only   in  the  area  AEBB'E'A'   where   the  Mohr- 
Coulomb  yield   criterion  applies; 

3)  The  active  earth  pressure  acts    on  plane  AA'; 

4)  Plane   strain   conditions   exist  with   respect   to 
depth; 

5)  The  piles   are   rigid  elements; 

6)  The   frictional   forces   on   surfaces   AEB  and  A'E'B' 
are   neglected  when  the  stress   distribution  in 
the   soil   AEBB'E'A'    is    considered. 

The   last  assumption  has   been  a  point   of   contro- 
versy  in  the  past,    because   it  seems   to  indicate  that 
the  stresses   acting  on  surfaces  AEB   and  A'E'B'    are 
principal  stresses.      This,    however,    is   not    true    (Ito 
and  Matsui,    1978);    Let  us   consider  two  small  elements 
within  the  plastic  region   (Figure  2),    element  I  in  the 
center  of    this   region,    and  element   II   on  the  edge.     The 
stresses   shown  in  element  I,    a&J  and   o^,    are  principal 
stresses,    and  the  Mohr's   circles  for  that  element  are 
shown  in  Figures   3a  and  3b  for  purely   cones ionless   and 
cohesive  soils,    respectively.   It   is   expected  that   the 
Mohr  circles  will  move  to  the   right  as   different   ele- 
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Element  I 


Figure  2      Stresses  Acting  On  Elements  in  Plastically 

Deforming  Ground  (After  Ito  and  Matsui,    1975), 
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Element  I   ^ 


Element  II 


Ca)  Cohesionless  Soil 


(b)  Frictionless   Soil 


Figure  3     Mohr  Circles  of  Elements  in  Plastically  Deforming 
Ground     (After    Ito  and   Matsui,    I  979). 
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ments  are  considered,  advancing  from  the  center  towards 
the  edge  of  area  GHH'G'  (i.e.,  in  Figure  3,  point  N 
will  move  along  the  failure  envelope  to  point  M).  How- 
ever, it  would  be  difficult  to  analyze  such  a  complex 
distribution  of  stresses.  For  simplicity,  every  point 
on  area  GHH'G',  except  planes  GH  and  G'H',  is  assumed 
to  be  under  the  same  stresses  as  element  I.  For  planes 
GH  and  G'H',  it  is  assumed  that  the  normal  stress  act- 
ing on  element  II,  o  -._,  is  equal  to  the  normal  stress 
acting  on  element  I,  o  _  (Figure  3).  It  is  clear  that 
the  assumption  of  a  uniform  stress  distribution  for 
AEBB'E'A'  does  not  consider  planes  EB  and  E'B'  to  be 
principal  planes.   The  stresses  acting  on  those  planes 
are  represented  by  points  M  and  M',  and  are  taken  into 
account  in  the  equilibrium  equations  written  for  area 
GHH'G'  (Equation  3-2  below). 

Derivation 

The  successive  steps  necessary  to  derive  the  force 
acting  on  the  pile  per  unit  depth  are  presented  herein. 
Expressions  of  this  force  are  given  for  c-«J>,  <J>,  and  c 
materials. 

First,  the  equilibrium  of  the  differential  element 
GHH'G'  surrounded  by  solid  lines  in  Figure  4a  is  con- 
sidered. Summation  of  forces  in  the  x  direction  gives: 
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Figure  4      Stresses  on  Elements  in  Plastically 
Deforming  Ground, 
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-Dda     -   a  dD  +  2dx   {a  tan   (^  +  -|)  +o  tan  9  +  c}   -  0(3-2) 


As   it   is   assumed  that    o  T,    is  approximately  equal  to 
a     ,    then   a  TT    corresponds   to  a  principal  stress    a 
from  assumption  6.     According  to  the  Mohr-Coulomb  yield 
criterion,    failure  occurs  when  the  stresses   satisfy   the 
equation: 


where 


o     -  oN     +  2c,  |N  .  (3-3) 

a         x   <J>  \  I    9 


N^  =  tan2    (£  +  |)  (3-4) 


From  geometry 


d(f) 

dx  =  -±— -  (3-5) 

tan(^  +  |) 


Substituting  Equations  (3-3)  and  (3-5)  into  (3-2): 


Ddo     =  dD    {N^tan^  N  -   l)a  +  c(2tan<j>  +  2N1/2+  N~1/2)} 
x  <J>  9x99 


(3-6) 


Integration  of    the  above  equation  gives    the  state   of 
stress    in  zone  EBB'E': 


(N1/2tan9+  N  -1) 
(C.D)      *                       *       -  c(2tan9+  2N1/2+  N   1/2) 
0      --i r75 * * 

X  NV    tan  9  +  N.-l 

9  9 


(3-7) 
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where  C.    is   a   constant    of    integration. 


Similarly,    the  equilibrium  of    the  differential 
element   shown  in  Figure   4b  will   give   the  state   of 
stress    in  zone  AEE'A'    (Figure   1).      Summing  all  forces 
in   the   x  direction: 


D„da     =■  2(o  tan    <f>  +  c)dx  (3-8) 

2     x  a 


Substituting  Equation    (3-3)    into  Equation    (3-8)   and 
integrating   it: 


2N  tan$  ,,- 

C„  exp    (— £ x)   -  c(2N  /    tan<jH-l) 

2  D  4> 

0     -  Z     M     -       . (3-9) 

x  N  ,  tan<b 


where  C_   is   a   constant   of    integration. 

Finally,    assumption  3  gives   the  normal  stress    on 
plane  AA'    (x  -   0): 


a  I      n  -   yzN"1  -   2cN'1/2  (3-10) 

x'x=0  <j»  4> 


where  z    is    the  depth  from  ground  surface  and   Y  Is   the 
unit  weight   of    the  soil.     Recognizing  that   the  plane 
AA'    is    the   x  =  0  plane,    Equation    (3-10)   can  be   con- 
sidered as  a   boundary   condition.     The   constant  C2  is 
obtained  by   introducing   this    boundary   condition  in 
Equation   (3-9): 
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C2   -   yz   tan<|i  +  c  (3-11) 

The   normal  stress   acting  on  plane  EE '   Is   found  by  using 
Equation   (3-9)  and   (3-11)  at   the   distance  x  -  AE. 

From  geometry, 


AE   -   {(Dx   -  D2)/2}   tan(ir/8  +   <fr/4)  (3-12) 


then, 


Klx=AE  "  NTi^  {<Y*tan«*c)-exp   (^  (3-13) 

N.    tan<J)  tan   (^  +-^))    -  c   (2N1/2tan4>  +   1)} 

9  Oh  <p 


The  constant  C.    is   obtained  by  considering  Equation 
(3-13)  as   a  boundary   condition  and  substituting  it   in 
Equation   (3-7). 


(N1/2tan<j*N  -1)    (N1/2tan<{rfN  -1) 
<ClV      *  *       -      Van/       {(Wn+fc) 


D  -D„  ... 

exp(-j> N  tan<|>  tan(Jf|))  -  c(2N^'   tan*  +  1)} 


1/2  -1/2 

+  c    (2tan<j>+  2N1/''  +  N//Z) 
?  9 


(3-14) 

The   lateral  force  per  unit   thickness   of   layer  acting  on 
plane  BB'   in  the   x-direction  is   obtained  by   substitut- 
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ing   the   above   value   in  Equation    (3-7): 


PBB'   "  VScW    "  A   [FT^  «Tf"«n*+c) 
1  9 


D1~D2  IT     A  1/2 

exp(-5— -N^tan*  tan(^-|))   -  c(2N7    tanf  +  1)} 


2tan<f>  +  2N1/2  +  N~1/2  2tan*  +  2N1/2  +  N~1/2 

+c J. f l  _  C*D i $ 

1/2                                J  1          1/2 

N   '    tan9+N-l  N^'    tan9  +  N^l 


where 


D,    (N1/2tan9  +  N  -1) 


(3-15) 


The   lateral  force  q   acting  on  the  pile  per  unit    thick- 
ness  of   layer  is    the  difference   between  the  forces   act- 
ing  on  planes   BB'   and  AA'    (Equations    (3-10)   and    (3- 
15)): 


1  VD2 

«  "  PBB'-D2(ax)x-0  "  A  C    [FTa^{eXP(-D7-%tan* 

9  * 


1/2  -1/2 

.  /0  2tan9  +  2NV      +  N  .  ' 

tan(|f|))   -   2N1/2tan9-l}H m * *— ] 


v  N/    tan*  +  N  -1 

9  9 


1/2  -1/2 

2  tan*+  2N7      +  N  .   '  .  /0 

-c{D. m x * 2D2N:1/2} 

1        N1/2tan9  +  N  -1  2    * 

9  9 
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4^2-  {A  exp(^— -  N^tan*  tan(|f|))  -  D2}  (3-16) 

<j»  2 


For  a  purely  cohesionless   soil,    the   cohesion  c  is 
zero  and  Equation   (3-5)   reduces    to   the   following 
expression  for   the   lateral  force  per  unit   thickness: 


-  D.    N^tan-trt-N  -1)  D-D. 

q  =  JE  {Di(   1)    *  *       exp(   li  H  tan*  t.n(^»-D2> 

*  2  2 

(3-17) 


For  a  purely  cohesive  soil  the  angle  of  internal 
friction  is  zero  and  several  of  the  above  equations  are 
modified,  Equation  (3-6)  becomes: 


,n   da 


The  state  of  stress  in  EBB'E'  is  obtained  by  integrat- 
ing this  equation: 

o  =  3c  logD  +  C3  (3-19) 

where  C,   is   a  constant   of   integration.     After  substi- 
tuting Equation    (3-3)   into  Equation    (3-8)   the  following 
expression  is   obtained: 


do  j 

— ^  =  jT-dx  (3-20) 

c         D2 


Integration  of  the  above  equation  gives  the  state  of 
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stress  in  AEE'A' 


Ox=^x+C4  (3-21) 


where  C,  is  the  constant  of  integration.  From  assump- 
tion 3  the  state  of  stress  in  plane  AA'  is: 

|ox'x-0  =  *~   2c  (3_22) 

If  Equation  (3-22)  is  considered  as  a  boundary  condi- 
tion, substitution  in  the  above  expression  gives  the 
following  for  C, : 


C4  -  Y2  -  2c  (3-23) 


Then  the  normal  stress   acting  on  plane  EE'   is: 

D1~D2 
laxlx={(DrD2)/2}    tan(w/8)   =  c(-^—  tan  £-2)+^    (3-24) 

Considering  Equation    (3-24)   as   a  boundary   condition,    C_ 
is   obtained   as: 


D-D 
C3   -  c   (-g tan|  -  3logD2  -  2)  +  yz        (3-25) 


Using  Equations    (3-25)   and   (3-19),    the  lateral  force, 
P„B,,    on   the   plane   BB'    is: 


V    -Dl{ax}D=Dl  <3"26> 
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Di     DrD2 

Dj    {c    (3log  ^  +  -g— =■  tan£  -   2)  +   72} 


Hence,    for   cohesive  soils,    the   lateral  force  q  per  unit 
thickness   of    layer  is: 


'-'bB'-V'x'x-O  (3"27) 


D    D-D 
c  {DjOlog  gi  +  ~-=-  tan|)  -  2^^)}  +  yz  (VV 


The  total  lateral  force  acting  on  a  stabilizing 
pile  due  to  the  plastically  deforming  layer  around  the 
pile,  F  ,  is  obtained  by  integrating  Equations  (3-16), 
(3-17)  and  (3-27),  along  the  depth  of  the  soil  layer 
for  c~4>,  <(>,  and  c  materials,  respectively. 

Parametric  Studies 

The  theory  of  plastic  deformation  involves 
numerous  parameters  that  describe  the  geometry  of  the 
problem  and  the  properties  of  the  soil.  As  seen  from 
the  final  expressions  3-16,  3-17  and  3-27,  these  param- 
eters are  the  unit  weight,  y»  the  depth,  2,  the  center 
to  center  distance  between  piles,  D. ,  the  clear  dis- 
tance between  piles,  D  ,   the  angle  of  internal  fric- 
tion, <)>,  and  the  cohesion,  c. 

The  lateral  force  per  unit  depth,  q,  is  a  linear 
function  of  parameters  Y  and  z.  Hence,  q  increases 
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linearly  as  y   or  z  increases.  The  effects  of  the 
remaining  parameters  are  more  difficult  to  detect.  The 
lateral  force,  F  ,  is  plotted  versus  the  ratio  ^/D^ 
for  different  values  of  c,  $,    and  pile  diameter,  b,  in 
Figures  5,  6,  and  7,  respectively.   In  general,  these 
results  indicate  that  the  lateral  force  decreases  as 
the  ratio  D2/Di  increases.   It  is  also  apparent  that 
the  force  increases  with  an  increase  in  4>  or  c.  This 
is  to  be  expected,  since  it  is  harder  for  a  stiffer 
soil  to  squeeze  between  the  piles.   Finally,  the 
lateral  force  increases  as  the  pile  diameter,  D.  -  D», 
increases  (Figure  7). 

Field  and  Laboratory  Measurements 

In  order  to  test  the  validity  of  the  theoretical 
approach  for  the  estimation  of  the  lateral  forces  on 
piles,  Ito  and  Matsui  (1975)  compared  observed  field 
values  to  calculated  ones.  The  observations  were  made 
in  three  different  locations  in  Japan,  at  Katamachi, 
Higashitono  and  Kaniyama.  Reinforced  concrete  piles 
about  1  ft  in  diameter  were  used  in  Katamachi  and  steel 
pipe  piles  of  about  the  same  diameter  were  used  at  the 
other  two  locations  for  the  stabilization  of 
landslides. 

The  field  lateral  forces  were  deduced  from  meas- 
ured strains  induced  in  the  piles  by  the  sliding  soil 
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masses.  The  theoretical  forces  were  estimated  using  c 
and  <J>  values  obtained  from  shear  and  standard  penetra- 
tion tests.  The  experimental  and  theoretical  results 
were  compared  and  several  conclusions  were  reached. 
Since  the  computed  and  measured  forces  were  found  to  be 
within  the  same  order  of  magnitude,  it  was  concluded 
that  the  theory  of  plastic  deformation  can  be  used  to 
predict  the  force  acting  on  stabilizing  piles.  When 
the  pile  head  was  restricted,  the  distribution  of  the 
force  was  trapezoidal,  which  agrees  with  the  theory. 
When  the  pile  head  was  allowed  to  deflect,  the  distri- 
bution of  the  force  was  triangular.   In  most  applica- 
tions, the  condition  of  restrained  pile  head  will  be 
used  and  the  force  on  the  pile  as  given  by  this  method 
will  give  reasonable  solutions.  As  will  be  seen  in  the 
remaining  of  this  work,  the  bending  moments  and  shear 
forces  acting  on  a  free  head  pile  are  much  larger  than 
those  acting  on  a  restrained  head  pile,  assuming  that 
the  forces  given  by  the  theory  of  plastic  deformation 
applies  to  both  piles. 

In  addition  to  the  field  tests,  Ito  and  Matsui 
(1982)  performed  a  series  of  model  tests  on  piles 
placed  in  sand  or  clay.  The  test  apparatus  consisted 
of  a  soil  container  with  the  model  piles  in  a  row  and  a 
lateral  loading  system.   The  soil  was  pushed  through 
the  piles  slowly  and  its  behavior  was  observed.   A 
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series  of  marks  were  placed  at  the  top  of  the  soil 
layer  to  indicate  the  shape  of  the  flow  lines  as  the 
soil  moved  between  the  piles.  It  was  observed  that  the 
disturbance  of  the  flow  lines  appeared  to  occur  mainly 
within  the  zone  of  plastic  deformation  that  was  origi- 
nally assumed.  This  validates  the  assumption  that  the 
soil  is  in  a  plastic  condition  around  the  piles. 

The  relationship  between  the  soil  displacement  and 
the  lateral  force  on  the  pile  was  also  investigated. 
It  was  found  that  the  lateral  force  increases  to  a 
yield  value,  and  then  reaches  an  ultimate  value  with 
Increasing  soil  displacement.   In  a  logarithmic  scale 
this  relationship  can  be  represented  by  a  bilinear 
curve  with  an  Inflection  point.  The  theoretical  lateral 
forces  given  by  Equations  3-16,  3-17,  and  3-27  are  very 
close  to  the  experimental  values  at  the  inflection 
point  (i.e.,  the  yield  value).  The  ultimate  lateral 
force  can  be  approximately  estimated  as  1.6  times  the 
theoretical  lateral  yield  force. 


SAFETY  FACTOR  OF  THE  STABILIZED  SLOPE 


The  stability  of  a  slope  can  be  investigated  by  a 
number  of  limit  equilibrium  methods,  including  the  log- 
arithmic spiral  method,  the  friction  circle  method,  and 
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the  method  of  slices.  Of  these,  the  friction  circle 
method  was  found  to  be  the  most  convenient  to  take  into 
account  the  force  exerted  on  the  slope  by  the  piles  in 
the  limit  equilibrium  analysis.  Parametric  studies 
were  performed  to  observe  the  change  in  the  potential 
failure  surface  with  the  addition  of  the  piles,  and  to 
obtain  the  relationships  between  the  safety  factor  of 
the  slope  and  the  location  and  spacing  of  the  piles. 

The  limit  equilibrium  calculations  are  based  on  an 
assumed  shape  of  the  rupture  surface.  The  safety  fac- 
tor, FS,  is  defined  as  the  ratio  of  the  shear  strength 
available  to  the  shear  strength  required  to  maintain 
the  slope  in  a  state  of  limit  equilibrium.  Assuming 
the  Mohr-Coulomb  failure  criterion,  the  factor  of 
safety  is  given  by: 


c     +  o  tan   <j> 

FS  =     a         n  -2.  (3-28) 

c     +   a  tan   <b 
r         n  Tr 


where   the  subscripts   denote  available  and  required 

quantities,    and   a    is   the  normal  force  acting  on  the 
n 

surface   of    rupture. 


As  an  aid  in  determining  FS,    the  factors   of   safety 

with   respect   to   cohesion,  F   ,    and  friction,   F  .,    can  be 

c  <p 

used: 


c  tan    $ 

F     -  —  and     F     = ■£  (3-29) 

c        c  *       tan    <b 

r  r 
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The  "true"  safety  factor,  FS,  is  obtained  when  F  and 


c 


F  are  equal  (Perloff  and  Baron,  1976): 


FS  =»  F  =  Fx  (3-30) 

c    9 


The  Friction  Circle  Method 

The  Friction  Circle  Method  (Taylor,  1937)  consid- 
ers the  stability  of  the  entire  sliding  mass  as  a  whole 
and  involves  taking  moments  with  respect  to  the  center 
of  a  circular  failure  arc.   By  investigating  all  poten- 
tial failure  surfaces  passing  through  the  slope,  the 
critical  failure  arc  and  safety  factor  are  evaluated. 
This  method  gives  reasonable  values  for  the  factor  of 
safety  when  compared  to  other  methods  of  slope  stabil- 
ity analysis   (Siegel,  1975).   When  all  the  normal 
stresses  acting  on  the  critical  surface  are  assumed  to 
be  concentrated  at  one  point,  a  lower  bound  value  for 
the  factor  of  safety  is  obtained.  When  the  normal 
stresses  are  assumed  to  act  at  the  two  end  points  of 
the  critical  surface,  an  upper  bound  factor  of  safety 
is  determined.   In  this  analysis,  the  normal  stresses 
are  assumed  to  be  concentrated  at  a  single  point  on  the 
failure  surface.  Therefore,   the  computed  factor  of 
safety  is  conservative. 
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The  model  used  by  Taylor  (1937)  for  the  analysis 
of  homogeneous  slopes  is  shown  in  Figure  8.   The  forces 
that  act  on  the  mass  are  the  weight,  W,  the  cohesion 
force  required  to  maintain  equilibrium,  Cr,  and  the 
resultant  of  the  normal  and  frictional  forces,  P. 

The  weight  W  acts  vertically  through  the  center  of 
gravity  of  the  mass  and  is  determined  by  multiplying 
the  unit  weight,  y.  of  the  soil  by  the  area  enclosed 
within  the  failure  surface  and  the  slope  boundaries. 

The  resultant  mobilized  cohesion  force  C   is  equal 
to  the  product  of  the  unit  cohesion  required  for 
equilibrium,  c  ,  by  the  length  of  the  chord  AB.   If  cf 
is  decomposed  into  its  components  parallel  and  perpen- 
dicular to  the  chord,  the  components  perpendicular  to 
the  chord  create  equal  and  opposite  moments  around 
point  0  (Figure  8).  Hence,  only  the  components  paral- 
lel to  the  chord  contribute  to  the  moment  equilibrium 
equation.  The  line  of  action  of  C  is  parallel  to  the 
chord  AB  while  its  moment  arm,  a,  is  obtained  by  equat- 
ing the  following  moments: 


c  .AB.a  =  c  .AB.R  (3-31) 

r         r 


Therefore, 


a  =  R.AB/AB  (3-32) 
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Figure  8     Forces  on  Unreinforced  Slope 
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The   force  P   is    the    resultant   of    the   normal  and 
frictlonal  forces   concentrated  at   one  point.      It   is 
almost   tangent   to  a  circle  of   radius  Rsin<|>  as   shown  in 
Figure  8. 

As  a  result   of   the  above  assumptions,    a  force 
polygon  can  be   constructed    (Figure   8).     A  detailed 
mathematical  solution  for   the   equilibrium  equations   of 
the  slope   is   given  by  Taylor   (1937).   The  following  two 
expressions   for   the  stability   number  were   derived  by 
Taylor   for  a   toe   failure  and  a   failure   below   the   toe, 
respectively: 


12  2 

c  -7?  esc  x  (ycsc  y  -  coty)  +  cotx  -  cotl 

a     -  ± - (3-33) 

F   tfl  2cotx  cotv  +  2 

c 


and 


c 


a 


12  2 

■sr  csc  x  (ycsc  y  -  coty)  +  cotx  -  coti  -  2r 


F   1&  2cotx  cotv  +  2 

c 


-(3-34) 


where 

x,y  =  angles   describing  the    <}rcircle   (Figure  8) 
i  ■  slope  angle 
v  ■  angle   defining  the  force  polygon  of 

the  slope    (Figure  8) 
H  ■  height   of    the  slope 
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n  =  ratio  of  the  distance  between  the  toe 
of  the  slope  and  point  A  (Figure  8) 
to  the  height  of  the  slope. 

The  safety  factor  with  respect  to  cohesion  on  any  sur- 
face defined  by  angles  x  and  y  can  be  obtained  using 
Equations  3-33  and  3-34.  The  true  safety  factor  of  any 
of  these  surfaces  is  obtained  through  successive  itera- 
tion until  F   is  equal  to  F  , .  The  critical  surface  is 
the  one  which  minimizes  the  factor  of  safety.   This 
minimum  value  is  the  safety  factor  of  the  slope. 

When  a  row  of  piles  is  inserted  in  the  slope,  the 
additional  resistance  provided  by  the  piles  changes 
both  the  safety  factor  and  the  critical  surface.  The 
new  slope  to  be  analyzed  is  shown  in  Figure  9.  The 
forces  applied  to  the  slope  are  identical  to  the  ones 
in  Figure  8,  with  the  exception  of  the  force  exerted  on 
the  slope  by  the  piles,  F  .  To  obtain  this  force  per 
unit  width  of  the  sliding  mass,  either  Equation  3-16, 
3-17  or  3-27  is  integrated  along  the  depth  of  the  pile, 

and  the  result  is  divided  by  the  center  to  center  dis- 

Ft 

tance,  D.  (i.e.,  F  ■  =— ). 
1      '   p   D  ' 

The   resistance  force,    F    ,    can  be  incorporated  in 

P 

the  force  polygon   (Figure  9)  resulting  in  two  new 
expressions  for  the  stability  number  for  a  toe  failure 
and  a  failure  below  the  toe,  respectively: 
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Figure    9     Forces  on  Slopes  Reinforced  with  Piles. 


E  - 


a    yH 


46 

12F 

p  fcos(CEO)  H 

■*-  { — j tt  cscx  cscy  sin*  +  OG} 

sin  v   2        '     T 


F  ^H   ,     2  r  cosx  ,     ,    ,.     /    \i 

c     o  esc  x  cxcy  sin®  {— : +  csc(u-v)  cosCx-u)} 

J  T  sin  v 


(3-35) 


where 


2 
E=l-2(cot  i  +  3coti  cotx  -  3coti  coty  +  3cotx  coty) 


and 


12F 

"a   _  7H" 


n 

{E  +  6n  -  6nsin<(cscx  cscy} ~^-  A 


(3-36) 


F  ^    c  2  •    r COSX         .     ,     ,     .  1 

c     ocsc  x  cscy  sin0 — : —  +  csc( u-v)cos (x-u ) \ 

lsinv 


where 

CEO  =  angle  F  forms  with  the  horizontal 

OG  =  moment  arm  of  F 

P 

u  =  angle  defining  the  force  polygon  of 
the  slope  (Figure  9) 


The  derivations  for  these  equations  are  given  in 
Appendix  A.  These  equations  can  be  used  in  exactly  the 
same  manner  as  Equations  3-33  and  3-34  to  obtain  the 
safety  factor  for  the  slope.   The  two  additional  param- 
eters are  the  force  F  and  the  angle  CEO,  which  defines 

the  direction  of  F  . 
P 


The  line  of  action  of  F  is  assumed  to  pass 

P 

through  the  centroid  of  the  force  diagram  given  by 


47 

Equations  3-16,  3-17  and  3-27.   This  assumption  is  rea- 
sonable since  each  individual  pile  is  preventing  move- 
ment essentially  by  cantilever  action.   The  direction 
of  F  is  assumed  to  be  parallel  to  the  line  tangent  to 
the  failure  surface  at  the  point  of  intersection  of 
that  surface  with  the  piles. 

The  magnitude  of  F  is  a  function  of  the  length  of 
the  pile  from  the  ground  surface  to  the  failure  sur- 
face, shown  by  distance  CE  in  Figure  9.   For  each  set 
of  angles,  x  and  y,  selected  to  investigate  a  particu- 
lar surface  in  the  slope,  the  distance  CE  changes  and, 
consequently,  the  magnitude  of  the  force  changes.   To 
take  this  into  account,  the  length  of  the  pile  above 
the  failure  surface  was  expressed  as  a  function  of  both 
the  ^circle  parameters,  angles  x  and  y,  and  the  loca- 
tion of  the  pile.  For  each  set  of  values,  x  and  y, 
chosen  to  determine  a  failure  surface,  a  new  pile 

length,  and  thus  a  new  F  ,  is  calculated.  This  force 

P 

is  then  used  in  Equations  3-35  or  3-36  to  determine  a 
new  stability  number.   The  calculations  necessary  to 
determine  the  length  CE  are  discussed  in  Appendix  B. 

Parametric  Studies 

A  computer  program  was   developed  to  perform  the 
two-dimensional  slope  stability  analysis  presented  in 
the  previous   section.      The   listing  and  user's   manual 
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for  this  program  are  given  in  Appendix  C.   It  can  be 
used  to  calculate  the  factor  of  safety  of  a  slope  made 
of  a  homogenous  material  and  reinforced  by  a  single  row 
of  piles.  The  soil  properties  are  defined  by  the  unit 
weight,  y.  cohesion,  c,  and  friction  angle,  <J>.   The 
program  was  used  to  perform  parametric  studies  to 
investigate  the  effect  of  the  location,  spacing,  and 
size  of  the  piles  on  the  stability  of  the  slope. 

When  piles  are  inserted  in  the  slope,  the  location 
of  the  critical  surface  changes  since  an  additional 
force,  F  ,  is  introduced  in  the  limit  equilibrium  equa- 
tions. This  is  illustrated  in  Figure  10  for  a  slope  of 
height  45  ft  and  angle  30  degrees  with  material  proper- 
ties c,  <J>,  and  y  equal  to  500  psf,  10  degrees  and  125 
pcf,  respectively.  The  original  factor  of  safety  of 
the  slope  (without  the  pile  reinforcement)  was  1.08, 
obtained  for  the  critical  surface  OAB.  After  insertion 
of  a  row  of  piles  with  diameter  ratio  D2/D.  of  0.6, 
placed  45  feet  upslope,  the  factor  of  safety  increased 
to  1.82  and  the  critical  surface  changed  to  O'AB'. 
Insertion  of  piles  70  feet  upslope  resulted  in  a  factor 
of  safety  of  1.64  and  the  critical  surface  0"AB". 

Without  the  piles,  any  surface  below  or  above  the 
critical  one  gives  a  higher  safety  factor.   With  the 
addition  of  the  piles,  every  possible  surface  in  the 
slope  is  reexamined.   To  obtain  a  factor  of  safety  for 
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each  surface,  the  weight  of  the  failure  wedge,  the  soil 

resistance,  and  the  reaction  from  the  piles,  F  ,  are 

balanced.  F  is  a  function  of  the  pile  length  from  the 

face  of  the  slope  to  the  surface  under  investigation. 

As  an  example,  when  the  surface  AEB,  which  was  the 

critical  surface  before  the  piles  were  inserted,  is 

examined,  the  distance  CE  is  used  to  calculate  the 

reaction  force  F  .  As  the  surface  becomes  shallower, 
P  ' 

this  distance  decreases  and  the  force  F  decreases. 

P 

As  F   decreases,  both  the  effect  of  F   on  the  safety 
P  P 

factor  and  the  rate  of  change  of  the  safety  factor 
decrease.  When  F  is  large,  it  has  a  predominant  effect 
on  the  safety  factor  whereas  the  other  resisting  forces 
are  negligible.  As  F  becomes  smaller,  the  effect  of 
these  forces  becomes  more  important.   Consequently,  the 
intersection  of  the  resulting  critical  surface  with  the 
piles  is  located  above  the  original  surface  (points  E' 
and  E"  in  Figure  10). 

Figure  11  shows  how  the  location  of  the  piles, 
indicated  by  the  distance  S  from  the  toe  of  the  slope, 
influences  the  factor  of  safety  of  the  same  slope  for  a 
given  ratio  D^/D..  For  each  value  of  S  the  safety  fac- 
tor was  computed  for  both  the  original  critical  surface 
(solid  curve  in  Figure  11)  and  the  new  critical  surface 
which  was  found  after  the  addition  of  the  piles  (dashed 
line  in  Figure  11). 
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These   results   indicate  that   the  safety  factor  com- 
puted on  the  original  surface   becomes  very  large  for  S 
between  45  and  75  feet.     For  this   range  of   values   of  S, 
the  force  exerted  by   the  piles   on  the  original  surface 
is   so   large   that   the   components   of   soil  strength  are 
not   required  for  equilibrium.      For   the  actual  critical 

surface,  the  reaction  force  F   is  smaller  and  the 

P 

actual  factor  of  safety  is  less  than  computed  for  the 

original  surface.   The  force  F  is  plotted  against  dis- 

P 

tance  S  for  both  the  original  and  the  new  surface  in 
Figure  12.  As  an  example,  for  S-45  feet,  the  force 
exerted  on  the  slope  (original  surface)  is  equal  to 
49,000  lb/ft,  and  the  factor  of  safety  2.10,  while  the 
actual  factor  of  safety  (critical  surface)  is  only 
1.82,  with  a  force  of  30,000  lb/ft. 

The  safety  factor  of  a  steep  slope  reinforced  with 
piles  reaches  a  maximum  when  the  piles  are  placed  very 
close  to  the  top  of  the  slope  (Figure  13).  The  criti- 
cal surfaces  of  a  steep  slope  remain  deep  and  the  fac- 
tor of  safety  keeps  increasing  as  S  increases  until  the 
piles  are  placed  very  close  to  the  top  of  the  slope 
(Figure  14).  Then,  the  surfaces  become  shallower,  the 
distance  CE'  is  very  small,  and  the  safety  factor 
starts  decreasing.  A  comparison  of  the  behavior  of 
steep  and  shallow  slopes  is  given  in  Figure  15.  The 

normalized  safety  factor,  FS/FS   ,  where  FS    is  the 

max         max 
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maximum  possible  safety  factor  for  a  constant  ratio 
D./D. ,  is  plotted  versus  the  normalized  distance,  S/L, 
where  L  is  the  horizontal  distance  from  the  toe  to  the 
top  to  the  slope.  The  figure  shows  that  the  piles  need 
to  be  placed  closer  to  the  top  of  a  steep  slope  than 
that  of  a  shallow  slope  for  a  maximum  safety  factor  to 
be  achieved. 

According  to  the  assumptions  made  by  Ito  and 
Matsui  (1975),  the  force  acting  on  the  slope  is  equal 
to  F  regardless  of  the  state  of  equilibrium  of  the 
slope.  Based  on  that  assumption,  the  stability  number 
can  be  expressed  as: 


Cfl/FcYH  -  f(F  )  (3-37) 


However,  an  overestimation  of  the  force  F  can  lead  to 

P 

unconservative  results  in  the  design  of  the  slope.  A 
more  practical  approach  for  design  is  to  introduce  the 

notion  of  a  mobilized  lateral  force,  F  ,  where: 

m 


F  -  F  /a  (3-38) 

m    p 


with  a  being  greater  than  1.0.   The  mobilized  force  is 
used  to  analyze  the  slope,  but  the  total  force  per  unit 
length  will  be  used  to  design  the  pile.   This  results 
in  a  conservative  design  for  both  the  piles  and  the 
slope. 
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It  is  proposed  herein  to  scale  the  force  F   by  the 
factor  of  safety  with  respect  to  cohesion  of  the  rein- 
forced slope  (i.e.  a  -  F  ).  The  resulting  stability 

c 

number  is: 


c  /F  1H   -  f(F  /F  )  (3-39) 

a  c        p  c 


This  equation  can  be  solved  by  iteration,  until  F  is 

equal  to  F.  (for  the  critical  surface,  F  -  F,  -  FS). 

<p  c  <p 

In  these  iterations,    the  force  F     is   divided  by   the 

P 

safety  factor,  FS  (i.e.,  F  =  F  /FS).   Physically,  this 

m     p 

implies  that  F  will  be  equal  to  F   for  a  slope  at  the 
m  p 

point  of  incipient  failure,  but  will  decrease  as  the 
original  degree  of  stability  of  the  slope  increases. 
The  critical  surface  falls  between  the  critical  surface 
obtained  without  the  piles,  and  the  one  obtained  with 
the  piles  providing  a  fully  mobilized  force.   The  curve 
relating  the  factor  of  safety  to  the  distance  S  (Figure 
16)  has  a  shape  similar  to  the  previous  curves.   How- 
ever, the  rate  of  increase  of  FS  is  less  than  for  a 
fully  mobilized  force,  and  the  peak  value  is  not  as 
sensitive  to  S  as  before.   For  the  example  in  Figure 
10,  piles  located  anywhere  between  45  and  65  feet  from 
the  toe  of  the  slope  provide  a  factor  of  safety  of 
about  1.55. 

The  previous  results  were  obtained  for  a  constant 
ratio  D-/D..   The  effect  of  this  ratio  on  the  factor  of 
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safety  is  illustrated  in  Figure  17  for  a  slope  of 
height  30  ft,  slope  angle  49  degrees,  and  a  distance  S 
of  7  ft.  As  expected  the  factor  of  safety  decreases 
significantly  as  D_/D.  increases. 

Figure  18  shows  the  relation  between  the  safety 
factor  and  the  distance  S  for  a  steep  slope  for  ratios 
D2/D,  of  0.538  (dashed  line),  and  0.6  (solid  line). 
The  force  per  unit  length  of  the  pile  given  by  Equa- 
tions (3-16)  and  (3-27)  decreases  with  an  increase  of 
the  ratio  D_/D. ,  regardless  of  the  position  of  the  pile 
upslope.  However,  to  find  the  force  F  ,  the  force  per 
unit  length  diagram  has  to  be  integrated  over  the  dis- 
tance CE  which  increases  as  S  increases.  Hence,  the 
difference  between  the  curves  becomes  larger  as  S 

increases  due  to  the  nonlinear  increase  of  F  with  S. 

P 

After  the  maximum  safety  factor  is  reached,  the  two 
curves  become  asymptotic  for  the  same  reason. 

The  variation  of  the  safety  factor  with  S  for  dif- 
ferent values  of  the  friction  angle,  $,  and  cohesion, 
c,  was  investigated  for  a  slope  of  height  18  ft,  angle 
30  degrees  and  a  ratio  Do/0]  °^  0.716.  The  shapes  of 
the  curves  (Figures  19  and  20)  are  the  same  as 
described  before.  Since  F  is  an  increasing  function 
of  $,  c  and  S,  its  influence  on  the  safety  factor  also 
increases  with  c  and  $. 
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Figure  I  7     Safety  Factor    versus  Ratio  D2/D|. 
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SUMMARY 


A  theoretical  approach  to  the  calculation  of  the 
lateral  forces  induced  on  piles  placed  in  moving  soils 
has  been  discussed.  Field  tests  indicate  the  the  pro- 
posed equations  can  be  used  to  predict  the  forces  on 
piles  used  in  the  stabilization  of  slopes.   Laboratory 
tests  also  validate  the  major  assumptions  behind  the 
theoretical  derivation.   The  Friction  Circle  Method  was 
adapted  to  incorporate  the  resistance  provided  by  the 
row  of  piles  into  the  slope  stability  analysis.   Com- 
puter programs  were  developed  to  compute  both  the 
lateral  force  and  the  factor  of  safety  of  the  slope. 
Parametric  studies  were  performed  to  investigate  the 
effects  of  the  geometric  and  material  parameters  on  the 
failure  surfaces  and  factors  of  safety  of  typical 
slopes. 
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CHAPTER  IV.   LATERALLY  LOADED  PILES 

In  designing  the  piles  to  resist  lateral  loads, 
the  deflection,  bending  moment  and  shear  force  profiles 
along  the  piles  are  required.   In  this  chapter,  these 
values  are  obtained  by  solving  the  differential  equa- 
tion (beam  equation)  governing  the  pile  displacements. 
A  closed  form  solution  of  this  equation  is  used  to 
analyze  the  pile  section  which  extends  above  the  criti- 
cal surface.  The  force  Intensity  on  that  section  is 
calculated  using  the  principle  of  plastic  deformation 
derived  in  Chapter  3  (Equations  3-16,  3-17,  and  3-27). 
A  finite  difference  method  is  proposed  to  analyze  the 
pile  section  which  is  embedded  below  the  critical  sur- 
face.  The  forces  acting  on  this  section  are  calculated 
using  the  subgrade  reaction  model.   Several  techniques 
to  obtain  the  modulus  of  subgrade  reaction  for  the  soil 
or  rock  below  the  critical  surface  are  discussed,  and 
recommendations  are  made  for  their  use  in  the  pile 
analysis  calculations. 
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THE  CONCEPT  OF  SUBGRADE  REACTION 

The  subgrade  reaction  model,  which  was  originally 
proposed  by  Winkler  in  1867,  characterizes  the  soil  as 
a  series  of  unconnected  linearly  elastic  springs 
(Hetenyi,  1946).   The  beam  (or  pile)  reaction  at  a 
point  is  simply  related  to  the  deflection  at  that 
point.   One  disadvantage  of  this  soil  model  is  the  lack 
of  continuity.   It  is  obvious  that  the  displacements  at 
a  point  are  influenced  by  stresses  and  forces  at  other 
points  within  the  soil.   However,  the  subgrade  reaction 
approach  provides  a  relatively  simple  and  efficient 
analysis  and  enables  factors  such  as  variation  of  soil 
stiffness  with  depth  and  layering  of  the  soil  profile 
to  be  taken  into  account.   It  is  used  here  to  determine 
the  displacements,  bending  moments  and  shear  forces  in 
a  pile  stabilizing  a  slope. 


GOVERNING  EQUATIONS 

A  beam  (or  pile)  supported  along  its  entire  length 
by  an  elastic  medium,  and  subjected  to  a  system  of  con- 
centrated forces  and  distributed  loads,  deflects  and 
creates  continuously  distributed  reaction  forces  In  the 
supporting  medium.   It  is  assumed  that  the  intensity  p 
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of  these  forces  is  proportional  to  the  deflection  y  of 
the  beam  (Figure  21a): 

P  -  Ky  (4-1) 

For  a  pile  of  width  b,  the  elastic  constant  K  is 
related  to  the  modulus  of  subgrade  reaction  K  by: 


K  =  b  Kg  (4-2) 


The   differential   equation   governing   the   beam  dis- 
placements,   y,    is   obtained  by   satisfying  the   equili- 
brium of   an  element    of    length   dz    (Figure  21b): 


.4 


HF-- 

dz 


EI  *-f  =    -Ky   +  q  (4-3) 


where:   E  =  Modulus  of  elasticity  of  the  pile 
I  =  Moment  of  inertia  of  the  pile 
q  =  Intensity  of  the  distributed  load 

The  geometry  of  the  pile  used  to  stabilize  a  slope 
is  shown  in  Figure  22.  The  following  notation  is 
adopted  in  this  figure  and  In  the  sequel: 


z  =  depth  from  the  ground  surface 
z  =  depth  from  the  sliding  surface 

q.  =  force/unit  length  acting  on  the  pile 
at  the  ground  surface  (z  =  -CE) 

q_  =  force/unit  length  acting  on  the  pile  at 
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Figure  2  I  (a)     Beam  on  Elastic  Foundation 
(After  Hetenyi,    I  946). 


M 


/^                 ^ 

S 

*— 

kydZ              *: 

"2          qd^ 

-< >- 

-* 

S  +  dS 


M  +  dM 


Figure   2  I  (b)    Cross  Section  of  a  Beam  on  Elastic  Foundation 
(After   Hetenyi,    I  946). 
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Figure   22     Stabilizing  Piles  Embedded  in  Bedrock, 
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the  failure  surface  (z  =  0) 
CE  ■  length  of  pile  from  ground  surface  to  failure 

surface 
ED  ■  length  of  pile  from  failure  surface  to  the  pile 

tip 
EB  =  length  of  pile  from  failure  surface  to  bedrock 

surface 
BD  =  length  of  pile  from  bedrock  surface  to  the 

pile  tip 

It    is    convenient   to   decompose   the   pile   equation 
into   two   equations    governing   the   pile   deflection  above 
and  below   the   failure   surface,    respectively: 


A, 

EI  — r=-  =  q(z)  (-CE    <  z    <  0)        (4-4a) 

dz 


A2 

EI  —^-  =   -K  y  (z    >  0)  (4-4b) 

dz 


where  y.  and  y„  are  the  pile  deflections  above  and 
below  the  sliding  surface  respectively,  and  q(z)  is  the 
applied  lateral  pressure  given  by: 


q2~ql 
q(z)  =  q2  +-£e-Lz  (*-5> 


where   q.    and  q~   are   obtained  according  to  the  methodol- 
ogy  developed   in  the  previous   chapter. 
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SOLUTION   TECHNIQUES 


Equation    (4-4a)   is   a  simple   beam  equation  and  a 
closed  form  solution  can  be   readily  obtained  by   direct 
integration.   A  finite   difference   scheme  was    chosen  to 
solve  Equation   (4-4b),    because   closed   form  solutions    to 
this   equation  are   difficult   to  develop  when  the  elastic 
constant    is    varying  with   depth. 

Closed  Form  Solution  of   Equation   (4-4a) 

The   solution  to  Equation   (4-4a)   is   obtained  by 
direct    integration: 


f  f      z5 

2  3  14  2 

yx   =  aQ  +  8lz  +  a^     +  a3z     +  j^^  z     +  y^^-gj     (4-6a) 


where 

f,    =  q. 


.  q2"qi 

2  CE 

and  an,    a.,    a„  >    and  a_   are   constants   of    integration. 

The  slope,    bending  moment,    and  shear  force   can  be 

obtained   by: 


dyx 

Slope   =  Ii- 

d2yx 

Moment   =  -EI  — ir-  (4-6b) 

dz 
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Shear  =  -EI  — r- 
dz 


where 


dYi  2    fl   3  ,    f2   4 

I  =  &1  +  2a2z  +  3a3z  +  ^-gj  z  +  ^Til  Z 


dz 


A1  f  f 

d  y,  fl   2     2   3 

— T  =  2a2  +  6a3z  +  2-EI  z  +  6H  Z 
dz 

and 


3 
d  yx        fx      f2   2 

T"  =  6a3  +  EI  2  +  2~EI  Z 
dz 

Sign  conventions  are  given  in  Figure  21b. 

Finite  Difference  Solution  of  Equation  (4-4b) 

A  finite  difference  scheme  is  used  to  solve  Equa- 
tion (4-4b).   The  embedded  length  below  the  failure 
surface  is  discretized  in  MT  equally  spaced  intervals 
of  length  X   (Figure  23).   The  interval  length  is: 

* ■ i  (W> 

and   the    location   of   any   node  m  is   defined  by: 

z   =  mX  (4-8) 


(4-6c) 
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Figure   23     Finite   Difference  Solution  for  a   Pile, 
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The  discretized  form  of  Equation  (4-4b)  is: 


A 

EI  (  &>m  =  b(Vm  ym  (4_9) 

,4m      o  m  m 
dz 


where   EI    is    the  pile   stiffness,    b  is    the  pile   diameter, 

(K„)      is    the   modulus    of    subgrade    reaction   of    the   soil 
o  m 

A 

at   node   m,   y      is    the   deflection  at   point   m,    and   ( — j-) 

m  dz* 

is    the   fourth   derivative    of    the    displacement   at   node   m. 


In  general,    the   modulus   of   subgrade    reaction 

varies  with   depth,    and   the   value    (Kc)      in  Equation    (4- 

b  m 

9)  is  different  for  each  node  m.   This  can  be  taken 
into  account  either  by  inputting  a  value  for  the 
modulus  at  each  node  or  by  introducing  an  empirical 
variation  with  depth.   Palmer  and  Thomson  (1948)  pro- 
posed the  following  expression  for  the  variation  of  K 
with  depth: 


KS  -  KSL  (E§>n  (4-10a) 


where  K   is  the  value  of  Kg  at  the  pile  tip  (z  =  ED) 
and  n  is  an  empirical  constant  equal  to  or  greater  than 
zero.   According  to  Poulos  and  Davis  (1980),  this  equa- 
tion Is  widely  used  to  predict  pile  behavior  using  the 
subgrade  reaction  approach.   It  is  proposed  herein  to 
modify  this  expression  to  account  for  soil  resistance 
at  the  point  z  =  0  (i.e.,  at  the  failure  surface). 
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Also,    two   different   equations  will   be  used   to   represent 
the   material  above   and  below   the  soil   rock    interface, 
respectively: 


h   ■   KS0  +   <KSL-   KS0>    <eI>"1  *    '  EB  (4-'°b) 


Ks "  So  +  (,Sl-  V  <ii2»n!     *  > EB        (4-10c> 


where  K       =   K     at   z   =  0 

KSL  "  KS  at  Z  "  EB 
1C  =  coefficient  of  subgrade  reaction  of  the 

material  below  z  =  EB 

K^  =  KR  at  z  =  EB 

KRL   =   KRatZ   =  ED 
n.    and  n„   are   empirical   indices    greater   than  or  equal 

to  zero.      Typical  values   for  K    ,   1C,    n.    and  n„  will   be 

presented  in   later  sections. 

These  expressions   for   the   soil  constant   can  be 
introduced  in   the   finite   difference   equations   to  pro- 
vide a  set   of   simultaneous   equations  with   the   nodal 
displacements   as   unknowns.      Finite   difference   approxi- 
mations  for   the  first,    second,    third  and  fourth   deriva- 
tives  of    the   displacement    below   the   critical   surface 
are: 


*[2.  ym+l  "  ym-l 

Cdz    V   =  2  A 
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( h     -  -2±i _E 5_L  (4-H) 

dz      m  .2 


,d3y2,  ym+2  -   2ynH-l  +  2ym-l  "  ym-2 


dz3     m  2X3 


and 


/y2,  ym-2  I   4ym-l  +   6ym  "  4ynH-l  +  ym+2 

\  4  ;m  A 

dz  A 


Substituting  equations  (4-11),  (4-10b)  and  (4-10c)  in 
the  finite  difference  equation  (4-9),  the  following 
expressions  are  obtained  for  each  node  m 

for  z  <  EB 

ym-2-4ym-l+6ym-4ym+l^nH-2  =  a  b[KS0+(KSL-KSO)(ET)ni]ym 

(4-12a) 

and  for  z  >  EB 

ym-2-4ym-l+6ym-4ym+l4ynH-2   =  a    ^^O^hL^RO^^^K 

(4-12b) 

where 


a   =  -^  (4-12c) 
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The  number  of  equations  that  the  finite  difference 
scheme  generates  is  equal  to  the  number  of  nodes,  MT  + 
1.   Four  additional  imaginary  points  (Figure  23)  are 
needed  for  the  complete  solution  of  equation  (4-9), 
which  brings  the  number  of  unknowns  to  MT  +  5.   In 
addition  to  the  MT  +  5  unknown  displacements,  four  con- 
stants of  integration,  an,  a,,  a-,  and  a,  are  unknowns. 
This  brings  the  total  number  of  unknowns  to  MT  +  9. 
The  total  number  of  equations  includes  the  MT  +  1  fin- 
ite difference  equations,  two  boundary  conditions  at 
the  pile  top,  two  boundary  conditions  at  the  pile  end 
and  four  continuity  equations  (at  the  failure  surface), 
which  brings  the  total  number  of  equations  to  MT  +  9, 
equal  to  the  number  of  unknowns.   The  boundary  condi- 
tions and  continuity  equations  are  discussed  in  the 
following  section. 

Boundary  Conditions 

To  ensure  continuity  of  the  pile  at  the  sliding 
surface,  the  following  continuity  relations  are  neces- 
sary: 


[y]z=o  =  [yi]z=o  =  [y2]z=o 


dyx  dy2 

t6]z=0   =    tdz_1z=0  =    [dz~ ]z=0 


a1  a1 

d  y,  d  y, 

[m]z=0  =  i-rh=o =  ittU     (4"13) 

dz  dz 
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tV]z=0  -  t-^]z=0  =  t-3-lz=0 
dz  dz 


where  9  is  the  slope  of  the  beam,  M  is  the  moment,  and 
V  is  the  shear  at  z  =  0. 

Using  the  above  equations,  the  constants  of 
integration  aQ,  a^  a2,  and  a3  are  expressed  as: 

a0  =  y0 


y1  -yx 


a,  = 


1     2X 


2   .  HZ^f^L  (4-U, 


6a,  = 


y2  -  2yj  +  2y_x  -  y_. 


3  2A3 


In  addition  to  the  continuity  equations,  boundary 
conditions  at  the  top  and  bottom  of  the  pile  have  to  be 
satisfied.   If  the  pile  is  driven  in  hard  rock  a  fixed 
end  condition  can  be  assumed;  In  this  case,  the  deflec- 
tion and  rotation  at  the  end  of  the  pile  (z  =  ED)  are 
prevented  (y  =  0,  6  =  0),  resulting  in: 


y   =0 

yMT 


and,    yMT+1  -  ym  -   1  -  0 


(4-15a) 
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If  the  pile  is  driven  in  soft  rock  the  moment  and 
shear  at  the  end  of  the  pile  can  be  assumed  to  be  zero. 
In  this  case,  the  two  boundary  conditions  become: 


MT+2 


yMT-2   "   2yMT+l   +   2yMT-l 


v  -   2v       +  V        ,    =  0  (4-15b) 

yMT+l  yMT        yMT-l  V 

At   the  top   of    the  pile    (z   =  -CE),    the   boundary 

conditions    depend  on   the   type   of  restraint.      One   of    the 

following  four   conditions   should  model  closely    the   res- 
traint  of    the   pile   in  the   field. 

Condition   1   -  Free   Head 

The   moment   and   the   shearing  force   at    the  pile  head 
(z   =  -CE)   are   zero: 


'M1Z~CE  •   "EI    t-rU-CE  "   ° 
dz 

and,  (4-16) 

*\ 

mz~CE--EI    l-T-'z-CE-0 
dz 


Substituting  equation  (4-6c)  into  the  above  equations 
gives : 
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2a2  -   6a3CE  +  ^±  CE*  -  ^|  CE3  =  0 

and,  (4-17) 


fl  f2        2 

6a3   "  EI   CE  +  2ET  CE     "  ° 


The  following  relations   are   obtained  for  the  moment   and 

shear   by    substituting   equation    (4-14)   in  equation    (4- 
17): 

(CE)y_2  +   (2A  -   2CE)y_1    -    (4A)yQ   +   (2  A  +  2CE)yx    -    (CE)y2 


3  2        3  12 

=   2XJ    [gg|  CEJ   -  ^eY  CE^]  (4-18a) 

and 


f  f 

-y_2  +  (2)y_r  (2)y!+  y2  ■  2x3tfl  CE  "  HI  ce2]  (4_18b> 


Condition  2  -  Unrotated  Head 

The  slope  and  the  shearing  force  at  the  pile  head 
(z  =  -CE)  are  zero. 


[6]z=-ce  =  i^r-WcE-0  (4"19a) 


[V]Z=-CE  =  "EI  [TT-]z=-CE  =  °  <4"19b> 

dz 


The  final  expression  for  the  shear  is  the  same  as 
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the  one  for  the  free  head  condition  (Equation  4-18b). 
The  expression  for  the  slope  is  obtained  by  substitut- 
ing equation  (4-6c)  into  equation  (4-19a): 


a^   -  2a2  CE  +  3a3  CE  -  ~:   CE  +  ^   |=  CE  -  0  (4-20a) 


Introducing  the   values   of   a.,    a„ ,    and  a~    given  in  equa- 
tion   (4-14),    the   following   relation  is   obtained: 


(-CE2)y_2+(2CE2-4  ACE-2    A2  )y_1+(8  XCE)yQ+(2  \2-4  ACE-2CE2)yi 


+   (CE2)y2  =   4  A3    [^gi  CE3  -  jjL.  CE4]  (4-20b) 


Condition  3  -  Hinged  Head 

In  the  hinged  head  condition,  the  deflection  and 
the  moment  at  the  pile  head  (z  =  -CE)  are  zero: 

^z=-CE  =  frll,-CB  =  °  (4"2U) 


[Mlz=-CE  =  "EI  hA-CE  =  °  (4"21b) 


The  final  expression  for  the  moment  is  the  same  as 
the  one  for  free  head  condition  of  equation  (4-18a). 
The  expression  for  the  deflection  is  obtained  by  sub- 
stituting equation  (4-6c)  into  equation  (4-2 la): 
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a0-aiCE  +  a2   CE2   -  a3   CE3  +  ^T  CE     "  T20~EI   CE     *   ° 

(4-22a) 


Introducing  aQ,    a^    a2,    and  a3   given  in  Equation    (4- 
14),    the   following   relation   is   obtained: 

(CE3)y_2  +   (6X2CE  +  6XCE2  -   2CE3)y_1  +   (12X3  -    12XCE2)y0 


+(6XCE2-6A2CE  +  2CE3)yi    -   (CE3)y2   -   12 X 


r L_  CE4  +  *—  CE5]  (4-22b) 

1    24EI  120   EI  J 


Condition  4  -   Fixed  Head 

The   deflection  and  the  slope  at   the   pile  head   (z   = 
-CE)  are  zero: 

fvl  =    [v    1  =0  (4-23a) 

lyJz=-CE        iylJz=-CE 


t^U-cE-^rWcE-0  (4-23b> 


The  final  expression  for  the  deflection  is  the 
same  as  the  one  for  the  hinged  head  condition  given  by 
equation  (4-22b),  and  the  expression  for  the  slope  is 
the  same  as  the  one  for  the  unrotated  head  condition 
given  by  equation  (4-20b). 
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COEFFICIENTS  OF  SUBGRADE  REACTION 


Coefficients  of  Subgrade  Reaction  for  Soils 

The  theory  of  subgrade  reaction  Is  used  to  model 
the  pile-soil  interaction  below  the  failure  surface. 
Estimates  of  the  coefficient  of  subgrade  reaction,  K  , 
are  necessary  to  evaluate  the  pile  displacements, 
moments  and  shear  forces.   In  this  section,  published 
criteria  and  data  on  K  are  reviewed  and  recommenda- 
tions are  made  for  their  use  in  design  applications. 

Two  fundamental  assumptions  are  made  in  this 
theory:   (1)  The  coefficient  of  subgrade  reaction,  K  , 
at  every  point  is  independent  of  the  contact  pressure, 
and  (2)  it  has  the  same  value  at  every  point  along  the 
contact  face.   Both  of  these  assumptions  are  approxima- 
tions of  the  true  conditions.   Loading  tests  performed 
on  actual  subgrades  show  that  the  settlement  does  not 
increase  linearly  with  increasing  pressure.  The  assump- 
tion of  linearity  is  usually  not  valid  for  pressures 
larger  than  about  half  the  ultimate  bearing  capacity  of 
the  soil.  This  limit  should  always  be  taken  into  con- 
sideration in  problems  involving  coefficients  of 
subgrade  reaction.   The  second  assumption  is  also  an 
approximation.  The  subgrade  reaction  does  not  have  the 
same  value  everywhere  on  the  surface  of  contact. 
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Depending  on   the   elastic   properties   of    the   subgrade, 
the  pressure   at   the    rim  is   either   greater   or  smaller 
than  at    the    center.      However,    the   errors    resulting   from 
this   assumption   can  be   neglected   in  practical  problems 
(Terzaghi,    1955). 

The   coefficient   of   subgrade    reaction  is    generally 
determined   by   one   of    the   following   three   methods:      (1) 
full  scale   lateral-loading  tests;    (2)   in  situ   tests 
such   as   plate    loading  tests,    pressuremeter  and   flat 
dilatometer   tests;    and,    (3)   empirical  correlations  with 
other  soil   properties. 

There   are   several   techniques    to  perform  lateral 
loading  tests    on  piles.      One   method   is    to  instrument 
the  pile  so  that   the  soil  pressures   and  the  pile 
deflections   are   measured   directly    (Matlock   and  Ripper- 
berger,    1958;    Bishop   and  Mason,    1954).      This   method 
provides   a   direct   evaluation   of  K   ,    but    it    is    time    con- 
suming and   relatively   expensive.      A  more   convenient 
procedure   is    to  measure   the   ground  line   deflection  and 
rotation  and   back-calculate  K     assuming  an  approximate 
distribution  with   depth    (Reese  and  Cox,    1969;    Welch, 
1972).      A  method  that  has    been  used  successfully   in  the 
past    is    the   placement   of    strain   gages    at   points   along 
the  pile.      The   strain   readings   are   converted   to  moment 
values   by   the  use   of    calibration   curves.      Then,    deflec- 
tions  are   computed      from        the   beam  equation. 
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The  use  of  plate  loading  tests  has  been  discussed 

by  Terzaghi  (1955).   The  main  problem  with  this 

approach  is  the  extrapolation  of  the  results.   Terzaghi 

assumed  that  the  coefficient  of  horizontal  reaction  for 

a  vertical  pile  surrounded  by  sand  is  a  function  of  the 

depth  below  the  surface,  z,  the  width  of  the  pile,  b, 

the  effective  unit  weight  of  the  sand,  y',  and  the 

relative  density  of  the  sand.   At  a  depth  z  below  the 

surface,  the  modulus  of  elasticity  of  sand,  E  ,  is 

s 

given  by: 


E  =  y'zA  (4-24) 

s 


where  A  is  a  dimensionless  coefficient  depending  only 
on  the  relative  density  of  the  sand  (Table  1).   Ter- 
zaghi proposed  the  following  relationship  between  K 

and  E  : 
s 


Ks  =  u33b  (4"25a) 


Substitution  of  equation  (4-25a)  into  equation  (4-24) 
gives: 


Ks  -  TTH  F  (4"25b> 


For  a  pile  embedded  in  stiff  clay,  Terzaghi  recom- 
mended a  constant  value  of  K  with  depth: 
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Table    1 
Values   of   dimensionless    coefficient  A,    to   calculate  kg 
(tons/ft    )   for  a  pile  embedded  in  moist   or   submerged 
sand    (Terzaghi,    1955). 


Relative   density   of    sand  Loose  Medium  Dense 


Range    of   values    of   A  100-300        300-1000        1000-2000 

Recommended  values   of  A  200  600  1500 


88 


KS  =  d~b  KS1  (4"26) 


where  K,,.    is    the   coefficient    of    vertical   subgrade    reac- 
tion  of   a  plate   of   width   equal  to  one   foot.      In  this 

equation   b   is   in  feet,    and   the   value    of  K     is   obtained 

3  — 

in  tons/ft    .      Recommended  values   of  K       are   given  in 

Table   2. 


Additional  suggestions  for  the  calculation  of  K 
include  the  equations  proposed  by  Chen  (1978),  which 
are  based  on  the  pressuremeter  modulus,  E  .   In  cohe- 
sionless  soil: 

E, 
Kg  =  3.3  f-  (4-27a) 

and  in  cohesive  soil: 

E, 
Kg  =  1.6  ^  (4-27b) 

2 
where  E,  is  the  pressuremeter  modulus  in  tons /ft  and 
d 

3 
K  is  given  in  tons /ft  . 


Another  relation  that  uses  the  pressuremeter 
modulus  was  derived  by  Yoshida  and  Yoshinaka  (1972): 


Ks  =  ^-Ed  (b)1/4  (4-28) 


Empirical  correlations  between  K  and  the  Young's 

Modulus,  E  ,  have  been  suggested  by  several  authors, 
s 

Vesic  (1961)  extended  the  work  of  Biot  (1937)  for  a 
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Table  2 

—  3 

Values  of  k  .  in  tons/ft  for  a  square  plate  (lxl  ft) 

resting  on  pre compressed  clay  (Terzaghi,  1955). 

Consistency  of  Clay         Stiff    Very  Stiff   Hard 

2 
Values  of  q   (tons/ft  ) 
u 

Range  for  k  .  (square  plates) 
Recomended  values  (square  plates) 


1-2 

2-4 

>4 

50-100 

100-200 

>200 

75 

150 

300 

90 

beam  placed  on  an  infinite  elastic  foundation  and  sub- 
ject to  a  concentrated  load.   His  results  indicated 
that  the  coefficient  of  subgrade  reaction  for  a  long 
beam  can  be  expressed  as: 


1  -  "s 


where  v     is  the  Poisson  ratio  and  E   is  the  Young's 
s  s 

modulus  of  the  subgrade.   Broms  (1964)  and  Francis 
(1964)  suggested  that  this  equation  can  be  applied  to 
piles,  and  they  used  it  to  estimate  the  modulus  K  for 
the  lateral  resistance  of  piles.   Francis,  however, 
considered  that  the  medium  extends  on  both  sides  of  the 
pile,  and  concluded  that  the  above  expression  should  be 
doubled: 


1  "  "s 


Equation  (4-29b)  can  be  used  for  either  sand  or 
clay  subgrades.   For  long  piles  in  soft  clay,  Gibson, 
in  an  unpublished  report,  suggested  that  a  soil  stiff- 
ness E  be  assessed  for  two  different  conditions, 

s 

instantaneous  and  long  term  loading.   Following  this 
suggestion,  E  may  be  calculated  as  either:   (1)  a 

5 

secant   modulus   measured  at   50%   of    the  ultimate 
undrained  triaxial   compression  test   for  instantaneous 


91 


loading  (u  should  be  taken  as  0.5  for  this  condition); 
s 

or  (2)  as  — —  for  the  long  term  condition  (u  ■ 

m  s 

v 

0.40),  where  m  is  the  modulus  of  volume  decrease 
v 

obtained  from  oedometer  or  drained  triaxial  tests.  For 
piles  in  sand,  any  of  the  equations  given  in  Table  3 
can  be  used  to  calculate  the  Young's  modulus.   These 
equations  have  been  used  to  predict  vertical  static 
compression  and  can  be  used  in  the  case  of  lateral 
loading  under  the  assumption  that  the  soil  is  isotropic 
and  homogeneous.   Correlations  between  E  and  the  cone 
resistance,  q  ,  need  some  additional  experimental 
verification  and  should  be  used  with  caution  (Jamiol- 
kowski  and  Garassino,  1977). 

In  general,  it  is  expected  that  the  coefficient  of 
subgrade  reaction  will  increase  with  depth  in  sand  and 
will  either  increase  or  remain  constant  with  depth  in 
clay.   Terzaghi  (1955)  recommends  values  of  nl  =  1  (to 
be  used  in  Equation  4-10)  in  sand  and  nl  =  0  in  clay. 
Davisson  and  Prakash  (1967)  suggested  that  nl  =  0.15  is 
a  more  realistic  value  for  clay  under  undrained  condi- 
tions.  Broms  (1974)  stated  that  the  coefficient  of 
subgrade  reaction  for  cohesive  soils  is  approximately 
proportional  to  the  unconfined  compressive  strength. 
As  the  unconfined  compressive  strength  of  normally  con- 
solidated clay  increases  linearly  with  depth,  the  coef- 
ficient of  subgrade  reaction  should  increase  in  a  simi- 
lar manner. 
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The  values  of  the  coefficient  of  subgrade  reaction 
can  be  obtained  by:   (1)  the  use  of  empirical  formulas 
(equations  4-25b  and  4-26);  (2)  the  use  of  in  situ 
parameters  (equations  4-27a  to  4-28);  and,  (3)  the  use 
of  empirical  correlations  (Equation  4-29b).   These 
equations  can  be  used  in  the  pile  analysis  computa- 
tions.  As  an  example,  assuming  that  the  material  below 
the  critical  surface  is  cohesionless,  either  equation 
(4-25b),  (4-27a)  or  (4-29b)  can  be  used  to  determine 
the  coefficient  of  subgrade  reaction.   When  equation 
(4-25b)  is  used,  the  constants  Kg0  and  KgL  needed  in 
equation  (4-12a)  can  be  obtained  by  equating  z  to  CE 
and  CB,  respectively  (Figure  22).  The  value  of  nl  for 
this  case  is  equal  to  1.   When  Equation  (4-27a)  is 

used,  a  value  of  (K_)   can  be  introduced  at  each  node  m 

'  am 

(m  =  0  to  m  =  MT,  as  shown  in  Figure  23).   When  Equa- 
tion (4-29b)  is  used,  the  values  of  Kg0,  KgL,  and  nl 
depend  on  the  equation  chosen  to  calculate  Eg  (Table 
3). 

Three  methods  to  calculate  the  modulus  of  subgrade 
reaction  in  soils  have  been  reviewed.  The  use  of  any 
of  these  methods  can  provide  satisfactory  results  since 
the  computed  pile  deflections,  bending  moments,  and 
other  quantities  will  not  be  very  sensitive  to  small 
changes  in  soil  modulus  values  (Matlock  and  Reese, 
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1960).   However,  the  use  of  pressuremeter  data  is 
preferable  since  it  eliminates  the  approximations 
required  to  interpret  plate  load  test  data  and  empiri- 
cal correlations.   In  addition,  the  pressuremeter  pro- 
vides modulus  values  in  the  horizontal  direction  which 
are  desirable  in  the  case  of  horizontally  loaded  piles. 

Coefficients  of  Subgrade  Reaction  for  Rocks 

A  soft    rock   can  be   treated  as   a  soil   in  the   sense 
that    it    can   be   given  a   value  IC   equivalent    to  K     in 
soil.      Geologists    overlook   this   value  and  prefer   to  use 
load-deformation   curves    to  predict   the   deformation 
modulus    or  in-situ  Young's  Modulus,   E    ,    of    the   rock. 
Therefore,    it    is   necessary   to  either  perform  field 
tests   or  use   correlations   between  1C    and  reported 
values    of    the  Young's   modulus. 

In  situ  methods   of    directly   obtaining  p-y   curves 
for   rock  are:      (1)    jacking  tests    (using   thin 
flatjacks);    (2)   dilatometer   tests;    and,    (3)   pressureme- 
ter tests.      The    jacking   test    is   similar  in  principle  to 
the  plate   load   test.      The   deformation   resulting  from 
the  application   of  a   load  of  known  intensity   to  a   given 
rock  area  is   measured  and  the   deformation  modulus   is 
calculated   by  using   the   theory   of   elasticity.   Perform- 
ing  jacking  tests   for  a  small   scale    job  is   not   practi- 
cal   (Wallace,    et   al.,    1970).      Data  from  such   tests   are 
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sometimes   available   for   rocks   frequently   encountered   in 
a   given   region.      Dilatometer   tests   have   several  advan- 
tages  over   the    jacking   tests:      (1)    the   dilatometer   can 
be   inserted  in  boreholes   of   conventional   diameter  for 
site   investigation;    (2)    the   tests   can   be   performed  at 
great   depths;    and,    (3)   large   radial  pressures   can  be 
applied    (Rocha,    1970).      Pressuremeter   tests   have  also 
been  used  successfully   in   the  past   for   the   calculation 
of    the   modulus    of   soft   rocks    (Dixon,    1970). 

A   load-deformation   curve   such   as    the   one   shown  in 
Figure   24   can   be   obtained   from  any   of    these   three 
tests.      The  modulus   of    this    curve   is   a  subgrade    reac- 
tion  modulus.      A   conservative   value   of    the   modulus    is 
obtained  by   taking   the  slope   of    the   curve   from  the   ori- 
gin to   the   point    of    total   deformation  and   maximum  load. 
Unfortunately,    these  values   are   rarely   reported  in  the 
literature,    but   are   merely   used   for   the    calculation   of 
the   deformation  modulus  E    . 

When  field  tests    cannot   be  performed,    K_   can  be 
back-calculated  using   the   values    reported   in   the 

literature   for   the   deformation  modulus,    E    .      For  exam- 

K 

pie,  if  the  rock  is  considered  elastic,  the  following 
relation  exists  for  a  modulus  obtained  with  a  jacking 
test: 


.    2 
y  =  p  b  C   I  U  )  (4-30) 

ER 
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Figure    24    Load  deformation  Curve  from  Plate  Jack  Test 
(After  Coon  and  Merritt,   1970). 
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where  b  Is  the  diameter  or  side  of  the  loading  plate 
and  p  is  the  pressure  applied  to  the  plate.   For  a  lxl 
ft  square  plate,  Equation  (4-30)  can  be  written  as: 


(4-31) 


y   (1  -  u2) 


If    the  in  situ  modulus,  E   ,    is   not   available,    it 
can   be   obtained   by  using   the    laboratory  Young's 
modulus,    E  _ -. ,    and  the  Rock.  Quality  Designation,   RQD. 
E   „    is    considered   to   be   the   characteristic  modulus   of 
the   intact   rock,    and  is    the   slope   of    the  stress-strain 
curve   at   50   percent    of    the  unconfined  strength.      Typi- 
cal values   of   E  __   are   given  in  Table   4.      It  has   been 

found   that   there   is    a   correlation   between  RQD   and  the 

ER 

ratio   of    the   moduli  =; (Figure   25  and  Table   5). 

t50 

These  empirical  relationships  are  based  upon  data  from 
twelve  construction  projects  (Coon  and  Merritt,  1970). 


SUMMARY 

A  method  to  analyze  piles  embedded  through  a 
failure  surface  into  either  bedrock  or  a  stiff  soil 
layer  has  been  developed.   The  portion  of  the  pile 
above  the  failure  surface  has  been  analyzed  using  a 
closed  form  solution.   The  portion  below  the  failure 
surface  has  been  analyzed  using  a  finite  difference 
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Table  4 
Modulus   of   elasticity   of   rocks. 


Rock  Modulus   of,Elasticity       Poisson's 

x   10      PSi  Ratio    » 


Granite  2.77-8.3  .15-. 24 

Gneiss  2.7-8.3  .11 

Quartzite         5.7-8.3 

Schist  .6-2.8  .08-. 20 

Sandstone  .69-11.1  .17 

Shale  1.4-4.9 

Mudstone  2.8-6.9 

Coal  1.4-2.8 

Limestone         1.4-11  .16-. 23 


103 


o 
m 


\ 

a: 
UJ 


< 
or 

CO 

3 
_) 
=) 
Q 
O 


1.0 


0.8  - 


0.6  - 


0.4  - 


0.2  - 


0 


•    Dworshak  Dam,  Granite  Gneiss,  Surface  Gages 
o  ii  "  "  "       ,  Buried 

o   Two  Forks  Damsite,  Gneiss 
■    Yellowtail  Dam,  Limestone 
a    Glen  Canyon  Dam,  Sandstone 

o   o  o 


y=0.023IX-l.32 


0 


20 


40 


60 


i 
80 


100 


ROCK  QUALITY  DESIGNATION,  ROD,  % 


Figure    25       Comparison  of  RQD  and  the  Modulus  Ratio  ER/Et50  (After  Coon 
and  Merritt,   I  970). 
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Table  5 
Correlation  between  def ormability  and  RQD 
(Coon  and  Merritt,  1970). 


Classification 

RQD 

Def  ormability  ■= 

Et50 

Very  Poor 

0-25 

>0.2 

Poor 

25-50 

>0.2 

Fair 

50-75 

0.2-0.5 

Good 

75-90 

0.5-0.8 

Excellent 

90-100 

0.8-1.0 
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solution  and  the  assumptions  of  the  theory  of  a  beam  on 
an  elastic  foundation.   Finally,  recommendations  for 
the  values  of  the  subgrade  reaction,  K  and  K  ,  as  well 

o        K 

as  their  variations  with  depth,  have  been  made. 

A  computer  program  which  uses  these  concepts  has 
been  developed  to  calculate  the  displacement,  bending 
moment  and  shear  force  along  the  length  of  the  pile. 
The  user's  manual,  a  listing  of  the  program,  and  an 
example  of  the  input  and  output  are  given  in  Appendix 
D. 
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CHAPTER  V.   SLOPE  STABILIZATION  AND  PILE  DESIGN 


Effective  stabilization  of  a  slope  with  piles 
requires  not  only  that  the  stability  of  the  slope  be 
assured,  but  also  that  the  piles  be  adequately  dimen- 
sioned.  In  this  section,  the  parameters  that  affect 
the  stability  of  the  slope  and  the  design  of  the  piles 
are  analyzed.  A  step-by-step  procedure  is  proposed  to 
select  these  parameters  and  achieve  an  efficient  sta- 
bilization scheme.  A  typical  application  is  given  for 
the  stabilization  of  the  shallow  slope  shown  in  Figure 
26.   The  factor  of  safety  of  the  slope  and  the  dis- 
placement, moment,  and  shear  profiles  along  the  piles 
are  calculated  by  using  the  methodologies  developed  in 
Chapters  3  and  4,  respectively.   Finally,  a  structural 
design  example  is  given  to  illustrate  how  various  fac- 
tors can  be  modified  to  arrive  at  an  optimum  pile 
design. 

The  slope  in  Figure  26  has  a  height  of  45  ft,  a 
slope  angle  of  30  degrees,  and  Is  made  of  a  homogeneous 
material.  The  friction  angle  of  the  soil,  <J>,  is  10 
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degrees,  the  cohesion,  c,  is  500  psf,  and  the  unit 
weight,  Y>  is  125  pcf.   The  critical  surface,  indicated 
by  a  dashed  line,  corresponds  to  a  minimum  safety  fac- 
tor of  1.08.   The  distance  from  the  ground  surface  to 
the  critical  surface,  CE,  is  19  ft.  A  factor  of  safety 
of  1.08  is  insufficient  and  the  slope  will  be  rein- 
forced with  a  row  of  piles.   The  following  steps  are 
proposed  to  achieve  an  efficient  design  of  the 
slope/pile  system. 

1)  A  parameter,  a,  which  represents  the  degree  of 
mobilization  of  the  force  F  ,  must  be  chosen.  As  indi- 
cated in  Chapter  3,  the  slope  stability  program  gives 

the  user  two  choices:   (1)  total  mobilization  of  F  ; 

P 

and  (2)  partial  mobilization  of  F  ,  in  which  case  a 

force  F  =  F  /a  is  used  to  represent  the  reaction  pro- 

m    p 

vided  by  the  piles.   It  is  proposed  herein  that  the 

ratio,  a,  between  the  total  force  and  the  mobilized 

force  be  taken  equal  to  the  safety  factor  of  the  slope, 

FS.  Physically,  this  implies  that  F   is  equal  to  the 

total  force,  F  ,  when  the  slope  is  in  a  state  of  limit 
P 

equilibrium,  but  it  decreases  as  the  degree  of  stabil- 
ity of  the  slope  Increases.  This  recommendation  will 
provide  a  conservative  assessment  of  the  stability  of 
the  reinforced  slope.   To  achieve  a  conservative  design 
of  the  piles,  the  notion  of  a  mobilized  force  is  not 
introduced  in  the  analysis  of  the  piles. 
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2)  The   horizontal   distance,    S,    between   the   pile 
row  and  the   toe   of    the  slope   may   be   dictated  by   site 
conditions   or  arbitrarily   chosen   (at    the  first  step   of 
the  design).      In  this   example,    S   is   taken  as   25   ft. 
With   the  piles   placed  at    that    location,    the   distance 
from  the   ground  surface   to  the   new   critical  surface, 
CE',    is   equal   to   16.5   ft. 

3)  The   slope   stability   program,    which   is    listed 
in  Appendix  C,    Is  used   to   find   the   factor   of   safety   of 
the   reinforced  slope,    as    a  function   of    the  pile   diame- 
ter,   b,    the   center-to-center   distance   between   the 
piles,    D. ,   and  the   location   of    the   pile   row  upslope, 
indicated  by   the  horizontal   distance   S  measured  from 
the   toe.      The   effect   of    both   the  pile   spacing  and   diam- 
eter  on  the  factor   of   safety   of    the  slope   can   be   con- 
veniently expressed  by  a  plot   of   the  factor  of   safety 
versus   the   ratio  D.  /b  for  a   given  value   of   S    (Figure 
27). 

4)  A  desirable  factor   of   safety  for  the  slope  is 
selected.      For   this   example,    it    is   assumed   that    the 
required  factor  of   safety   is    1.30,   which  is   a  value 
commonly  used   in   the  stability   of    slopes. 

5)  A  ratio  D  /b  can  be  selected  from  Figure  27. 
In  this  example,  the  ratio  should  be  less  than  2.8  to 
satisfy   the    requirement    of   a  minimum  factor   of   safety 
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Figure  27     Safety  Factor  versus  Ratio  D|  /b, 


Ill 


of  1.30  for  the  slope.   A  conservative  value  of  2.5  is 
chosen  here.  Note  that  large  ratios  D./b  should  be 
avoided  because  the  assumption  of  a  plastic  state 
around  the  piles  is  not  fulfilled  for  excessive  spacing 
between  the  piles. 

6)  In  the  initial  design  stage,  a  pile  diameter 
is  arbitrarily  chosen.   In  the  present  example,  the 
piles  are  assumed  to  be  2  ft  in  diameter  with  a 
center-to-center  distance  of  5  ft  (to  satisfy  a  D^/b 
ratio  of  2.5).   Assuming  that  a  reinforced  concrete 
pile  is  used,  the  pile  stiffness  can  be  expressed  as 
the  product  of  the  Young's  Modulus  of  concrete  and  the 
gross  moment  of  inertia  of  the  cross  section.   For  a 

concrete  with  a  compressive  strength  of  4000  psi,  the 

8      2 
pile  stiffness  is  4  x  10   lb-ft  . 

7)  The  pile  stability  program  listed  in  Appendix 
D  is  used  to  find  the  displacement,  bending  moment,  and 
shear  force  along  the  length  of  the  pile.   The  pile  is 
analyzed  as  a  beam  embedded  in  an  elastic  foundation. 
The  magnitude  of  the  force  per  unit  length  acting  on 
the  pile  section  above  the  critical  surface  is  obtained 
from  Equations  3-16,  3-17  or  3-27.  The  force  per  unit 
length  acting  on  the  pile  section  below  the  critical 
surface  is  a  function  of  both  the  pile  stiffness  and 
the  nature  of  the  foundation.   In  this  example,  the 
foundation  consists  of  a  clay  layer,  which  extends  10 
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ft  below  the  critical  surface  and  is  underlain  by 

bedrock.  The  coefficient  of  subgrade  reaction  of  the 

4 
clay  is  assumed  to  be  5  x  10  pcf  (Equation  4-26  and 

Table  2).   The  bedrock  is  a  soft  shale  with  a  coeffi- 

o 

cient  of  subgrade  reaction  of  4  x  10  pcf  (Equation  4- 
31  and  Table  4). 


The  displacement,  bending  moment,  and  shear  pro- 
files corresponding  to  the  selected  parameters  are 
given  in  Figures  28,  29  and  30  for  the  four  possible 
boundary  conditions  at  the  pile  top  discussed  in 
Chapter  4  (BC  =  1,2,3,4):   (1)  free  head;  (2)  unrotated 
head;  (3)  hinged  head;  and  (4)  fixed  head.   Figure  29 
reveals  that  the  hinged  head  condition  results  in  the 
smallest  bending  moment  in  the  pile,  followed  in  order 
by  conditions  (4),  (2),  and  (1).   Based  on  this,  a  res- 
trained pile  head  (hinged  or  fixed)  is  recommended.   In 
addition,  experimental  results  (Ito  and  Matsui,  1975) 
indicate  that  the  lateral  load  acting  on     piles  due 
to  plastically  deforming  ground  can  be  best  estimated 
by  the  theory  of  plastic  deformation  under  the  condi- 
tion of  a  restrained  pile  top.   Therefore,  the  free 
head  condition  should  be  avoided  in  order  to  both  limit 
the  moment  and  shear  on  the  pile  and  closely  estimate 
the  force  acting  on  the  piles.  A  restrained  head  con- 
dition can  be  obtained  by  connecting  the  pile  heads 
with  a  beam  which  is  fixed  by  tension  anchors.   The 
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unrotated  head  condition  can  be  obtained  by  simply  con- 
necting the  pile  heads  with  a  beam.   The  constraint 
used  for  the  pile  analysis  should  approximate  the  con- 
ditions in  the  field  as  closely  as  possible.   In  this 
example,  the  fixed  head  condition  is  assumed. 

8)  The  structural  analysis  of  the  pile  can  now  be 
performed.   The  maximum  displacement,  moment,  and 
shear,  acting  on  the  pile  are  the  three  factors  that 
should  be  considered  to  assure  that  the  design  is  ade- 
quate. 

First,  the  maximum  displacement  is  checked.   Shan- 
non and  Wilson  (1964)  used  a  limiting  displacement  of  2 
in.  to  design  the  piles  in  the  Seattle  Freeway  stabili- 
zation, and  this  value  is  assumed  here  as  the  maximum 
tolerable  deflection.   Choosing  an  appropriate  pile 
diameter  will  satisfy  this  restriction.   The  maximum 
displacement  for  the  fixed  head  condition  in  Figure  28 
is  0.24  in.,  which  is  well  within  the  tolerable  limit. 

Second,  the  steel  reinforcement  required  for  the 
pile  to  resist  the  maximum  bending  moment  is  designed. 
The  bending  moment  diagram  shown  in  Figure  29  has  a 
maximum  value  of  317  ft-kips  (boundary  condition  4). 
In  accordance  with  the  ACI  code,  the  moment  should  be 
increased  by  a  factor  of  1.7  to  a  design  value  of  539 
ft-kips.  The  CRSI  handbook  (1978)  is  consulted  to 
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obtain  the  appropriate  steel  design.   In  this  example, 
the  maximum  moment  can  be  resisted  using  12  No.  10 
rebars  of  grade  60  steel. 

Third,  the  maximum  shear  is  checked.   When  the 
pile  is  embedded  in  rock,  high  shear  forces  will  occur. 
The  shear  provided  by  the  concrete,  V  ,  is  considered 
to  be  the  product  of  an  average  shear  stress   (usually 
2\|f'c'  where  f'c  is  the  strength  of  the  concrete)  by 
the  effective  cross  sectional  area  of  the  member 
(McCormac,  1978).   In  this  example, 


V  =  2\|4000  x  452  =  5.7  x  1<T  lbs 


c 


The   maximum  shear   that   the   member  must  withstand  is    1.7 
times    the   maximum  applied  shear   obtained  from  Figure 
30,    and   Is   found   to   be   1.7  x   (7.3   x   10      lbs)   =    1.2   x 
10     lbs.      Since  this   exceeds    the  strength  provided  by 
the   concrete  alone,    stirrups   must   be  provided. 

As   an  alternative   design,    a  pile   of   3   ft   in  diame- 
ter  could  be  used.      In  this    case,    the   center-to-center 

spacing  is   increased  to  7.5  ft,    and  the   resulting 

9  2 

stiffness    is   2  x   10     lb-ft    .      The  maximum  displacement 

on  the   pile   is   0.1    inch    (Figure  31).      The   maximum  bend- 
ing moment    is   550   ft-kips    (Figure   32).      The  pile   can  be 
designed  with   12   No.    10   bars    to   resist   a   factored 
moment   of   935   ft-kips.      The  shear  strength   provided  by 
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the  concrete  is  1.29  x  10   lbs  and  the  factored  maximum 
shear  that  the  member  must  withstand  is  1.7  x  (1.1  x 
10  )  -  1.87  x  10   (Figure  33).   Stirrups  would  again  be 
provided. 

9)  An  optimum  design  can  be  obtained  by  minimiz- 
ing the  cost  of  materials  and  construction  for  dif- 
ferent configurations  of  the  slope/pile  system.  The 
two  parameters  that  can  be  varied  are  the  pile  diame- 
ter, b,  and  the  distance  of  the  pile  row  from  the  toe 
of  the  slope,  S.   For  this  example,  increasing  the  pile 
diameter  from  2  ft  to  3  ft  increases  the  required  spac- 
ing between  the  piles  from  5  ft  to  7.5  ft.   Hence,  for 
a  landslide  about  1000  ft  wide  ,  200  piles  of  2  ft  in 
diameter  or  134  piles  of  3  ft  in  diameter  will  assure  a 
factor  of  safety  of  1.30  for  the  slope.   Changing  the 
distance  between  the  pile  row  and  the  toe  of  the  slope 
has  also  an  effect  on  the  number  of  piles.   To  illus- 
trate this,  the  factor  of  safety  is  plotted  against  the 
ratio  D./b  for  a  distance  S  =  40  ft  upslope  (broken 
line  in  Figure  27).   In  this  case,  to  achieve  a  safety 
factor  of  1.30,  the  ratio  D./b  must  be  equal  to  3.5. 
Therefore,  only  143  piles  of  2  ft  in  diameter  or  100 
piles  of  3  ft  in  diameter  are  required.  However,  plac- 
ing the  piles  15  feet  further  upslope  increases  both 


*  A  slide  of  such  magnitude,  but  with  different  slope 
configuration  and  material,  occurred  during 
construction  of  1-471  in  Ohio. 
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the  pile  length  required  to  penetrate  the  critical  sur- 
face and  the  loads  that  act  on  the  pile,  and  conse- 
quently, the  required  reinforcement.   Hence,  this 
alternative  could  prove  to  be  less  economical  than  the 
first  one.   Other  considerations  in  the  selection  of 
the  final  design  would  include  the  degree  of  difficulty 
in  the  installation  of  the  pile  row  at  each  location 
and  labor  costs. 

10)   After  dimensioning  the  diameter  of  the  pile, 
the  length  of  the  pile  is  determined.   It  is  suggested 
that  the  pile  be  embedded  to  a  sufficient  depth  so  that 
the  bending  moment  and  shear  force  approach  zero.   To 
find  the  appropriate  depth,  a  pile  of  infinite  length 
is  analyzed  and  the  point  at  which  these  values 
approach  zero  is  located.   Embedding  the  pile  deeper 
than  this  point  will  not  increase  its  stability.  Con- 
sidering the  case  of  a  3  ft  diameter  pile,  the  bending 
moment  and  shear  force  are  zero  at  a  distance  18  ft 
below  the  critical  surface  (Figure  32  and  33).  There- 
fore, the  pile  will  penetrate  about  26.5  ft  of  soil 
(CE'  +  10  ft)  and  will  be  embedded  8  ft  in  the  soft 
rock.   Although  the  embedment  is  relatively  important 
for  this  particular  case,  it  has  to  be  recognized  that 
the  required  embedment  length  decreases  considerably  as 
the  stiffness  of  the  foundation  material  increases. 
The  piles  must  always  penetrate  the  critical  surface  of 
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the  unreinforced  slope.  Otherwise,  the  slope  would 
fail  on  its  original  failure  surface  and  the  piles 
would  not  contribute  to  its  stability. 

SUMMARY 

A  step-by-step  procedure  has  been  proposed  to 
analyze  the  stabilization  of  a  slope  with  a  row  of 
piles: 

1)  A  parameter  a  is  used  to  determine  the 
degree  of  mobilization  of  F  .   It 

is  proposed  herein  to  scale  the  force 
acting  on  the  piles  according  to  the  degree 
of  stability  of  the  slope  (a  ■  FS). 

2)  The  horizontal  distance,  S,  which  indicates 
the  location  of  the  piles  upslope  is 
selected. 

3)  A  graph  of  the  factor  of  safety  of  the  slope 
versus  the  ratio  D./b  is  constructed  using 
the  slope  stability  program  listed  in 
Appendix  C. 

4)  A  minimum  factor  of  safety  is  speci- 
fied. 


5)  The  ratio  D./b  which  satisfies  the  required 
factor  of  safety  is  selected. 
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6)  The  pile  diameter  b  is  arbitrarily  chosen, 
while  the  distance  D.  is  determined  by  the 
ratio  D./b. 

7)  A  boundary  condition  for  the  pile  head  is 
chosen.   This  condition  will  have  to  be 
fulfilled  during  construction. 

8)  The  displacement,  bending  moment,  and  shear 
force  along  the  pile  are  calculated  using 

the  pile  stability  program  listed  in  Appendix 
D.  Then,  the  pile  is  designed  structurally. 

9)  If  desired,  an  optimum  solution  can  be  obtained 
by  changing  the  pile  diameter  and/or  the  pile 
row  location  (i.e.,  the  distance  S) 

10)  The  depth  of  embedment  of  the  pile  is 
determined. 

This  systematic  approach,  used  in  conjunction  with  the 
programs  listed  in  Appendix  C  and  D,  provides  the  means 
for  an  efficient  selection  of  the  appropriate  pile 
dimensions  and  reinforcement  required  for  stabilization 
of  a  slope. 
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CHAPTER  VI.   SUMMARY  AND  CONCLUSIONS 

This  study  is  part  of  a  project  undertaken  at  Pur- 
due University  to  develop  a  methodology  for  the  design 
and  analysis  of  slopes  stabilized  with  piles.   Dif- 
ferent aspects  of  this  problem  have  been  considered, 
the  most  important  of  which  are:  (1)  the  calculation  of 
the  force  exerted  on  the  piles  by  the  slope;  (2)  the 
effect  of  a  row  of  piles  on  the  stability  of  the  slope; 
and  (3)  simultaneous  slope  stability  analysis  and  pile 
design  (dimensioning)  to  meet  minimum  safety  require- 
ments for  both  the  slope  and  the  piles.  The  following 
steps  have  been  taken: 

-Available  techniques  for  calculating  the 
pressure  acting  on  passive  piles  have  been 
reviewed.   The  "plastic  state"  method  (Ito 
and  Matsui,  1975)  was  selected  and  im- 
plemented in  a  computer  program.   This  pro- 
vides an  estimate  of  the  pressure  acting 
on  the  pile  regardless  of  the  state  of 
equilibrium  of  the  slope. 
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-A  computer  program  has  been  developed  to 
calculate  the  safety  factor  of  a 
slope  reinforced  with  a  single  row  of  piles. 
The  (jrcircle  slope  stability  method  has 
been  extended  to  take  into  account  the  force 
exerted  by  the  piles  on  the  slope.   It 
is  recommended  that  the  force  exerted  by 
the  piles  be  assumed  equal  to  the  maximum 
value,  calculated  using  the  plastic  state 
assumption,  divided  by  the  factor  of 
safety  of  the  slope. 

-Parametric  studies  have  been  performed  to 
obtain  relationships  between  the  safety 
factor  and  parameters  such  as  pile  diameter, 
center-to-center  distance,  and  location  of 
the  pile  row. 

-A  computer  program  has  been  developed  to 
compute  the  displacement,  bending  moment,  and 
shear  force  at  each  point  along 
the  pile.  The  pile  is  analyzed  in  two 
sections.   Since  the  pressure  acting  on 
the  section  above  the  critical  surface  is 
known,  this  section  is  solved  using  a 
closed-form  solution.   The  section  below 
the  critical  surface  is  analyzed  using 
the  finite  difference  method.   It  is  assumed 
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that    the  pressure  acting  at   any   point   along 
the  pile  is   a  function  of   the  deflection 
and   modulus    of    subgrade   reaction  at 
that  point.     Experimental  methods  which 
can  be  used  to   determine  the  modulus   of 
subgrade   reaction  have  been  discussed. 
Published  empirical   values  have   been   re- 
viewed. 

-A  step-by-step  procedure  has   been  proposed  for 
the   design   of   both   the   slope  and   the   piles. 
When  a   desired  safety   factor   for   the 
slope   and  a    location   of    the  pile   row 
are   chosen,    the  pile   diameter,    center-to- 
center   distance,    and  the   reinforcement   can  be 
determined  so  that   both   the  stability   of    the 
slope  and  the  integrity   of    the  piles   are 
insured.     A  typical  example,    including  the 
structural   design  of    the  piles,    has   been  given. 

The  most  important   conclusions   of   this   study  are 
as   follows: 

-After  installation  of   the  piles,    the  poten- 
tial failure  surface  is   expected  to  become 
shallower. 
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-As    the   distance  S,    which   indicates    the 
pile   location  upslope,    increases,    the 
factor  of   safety   changes  with  a   rate 
which  is  a  function  of   the  ratio  D2/Di 
(clear   distance/center-to-center   dis- 
tance  between  the  piles).      As   this    ratio 
decreases,    the    rate   of    change  in   the 
safety   factor   increases. 

-The  factor   of   safety   increases   more    ra- 
pidly  in  strong  soils    than  in  weak   ones. 
For   conventional   values    of    the   factor 
of   safety,    the  piles    can  be  placed  closer 
to   the   toe   for  a  strong   soil   than  for 
a  weak.   soil. 

-The  factor   of   safety   of   the  slope   de- 
creases  as    the   ratio   of    the   center-to- 
center  spacing  to  the   diameter  of   the 
piles,   D   /b,    increases. 

-For  a   given  ratio  D./b,    there   is   a 
range   of   values   of   S  which  will   satisfy 
a   required  factor   of   safety.      Generally, 
the  piles   must   be   located  closer  to  the 
top   of    the  steeper  slopes    than   of    the 
shallower  ones. 
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-The  pile  top  should  be  restrained  (fixed 
or  hinged  end)  to  minimize  the  bending 
moments  and  shear  forces  applied  to  the 
piles. 

-A  satisfactory  design  of  the  slope/ 
pile  system  can  be  achieved  to  satisfy 
the  stability  of  the  slope  and  to  in- 
sure the  structural  integrity  of  the 
piles. 
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APPENDIX  A 


Determination  of  a  Stability  Number  for  the  Reinforced  Slope 

Figure  9  in  Chapter  III  is  similar  to  the  figure 
used  by  Taylor  (1937)  for  the  derivation  of  a  stability 
number  for  an  unreinforced  slope.  Assuming  that  the 
force  F  is  known  in  magnitude  and  direction,  and  that 
its  presence  does  not  influence  the  point  of  action  of 
P,  a  new  stability  number  can  be  determined  for  the 
reinforced  slope. 

The  following  notation  is  used  in  the  derivation 
(Figure  9): 

W  =  weight  of  the  sliding  mass 
C  =  cohesion  required  for  equilibrium 

P  =  resultant  of  the  normal  and  frictional  forces 
F  =  pile  reaction  on  the  slope 

d  =  moment  arm  of  W 

a  ■  moment  arm  of  C 

r 

OG  =  moment  arm  of  F 

P 

H  ■  height  of  the  slope 
i  ■  slope  angle 
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x  =  counterclockwise  angle  from  the  horizontal 
to  the  chord 
2y  =  angle  AOB 
R  =  radius  of  the  circle 

Taylor  (1937)  derived  the  following  expressions  to 
arrive  at  a  stability  number.   The  cohesion  force 
required  for  stability  is: 


c  „ 

C  =c  AB  =  =^- AB  =  ^  c  R'sin(y)     (A-l) 
r    r      r       r   r 
c       c 


where 


c     =  unit   cohesion   required  for  equilibrium 

c     =  unit   cohesion  available 
a 

AB  =  length   of   the   chord 

F     =  safety   factor  with   respect    to   cohesion 


The   length   00 '    is    the  moment   arm  of   C     and  can   be 
expressed  as: 

00'   =  R  sin  (<())   csc(u-v)  (A-2) 

The  weight  W   is: 

W  =   fR2y-7R2    sin(y)   cos(y)+^-  {cot(x)-cot(i)}  (A-3) 

The   radius   of    the    circle   is: 


R  =  y  csc(x)   csc(y)  (A-4) 
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The  moment  of  the  weight  with  respect  to  the  center  of 
rotation,  0,  is: 

3 
Wd  =  -j|-{l-2cot2(i)+3cot(i)cot(x)-3cot(i)cot(y)   (A-5) 

+  3cot(x)cot(y)} 

For  the  sake  of  simplicity,  the  expression  inside  the 
parenthesis  will  be  replaced  by  E: 


Wd  "  U~  E  (A"6) 


Using  the  equations  derived  above  and  two  equa- 
tions of  equilibrium  (derived  by  summing  moments  around 
point  0  and  summing  forces  in  the  x  direction),  an 
expression  for  the  stability  number  is  derived: 


12  2 

c  ■=■  esc   (x)    {y  esc   (y)-cot(y)}+cot(x)-cot(i) 

(A-7) 


F   tfl  2  cot(x)   cot(v)  +  2 

where  the  angle  v  is   obtained  from  the  following  rela- 
tions: 


u 

sin(u-v)   -  ttt  sin(u)   csc(x)   csc(y)   sin   (<|>)      (A-8) 
with 


H  2 

cot(u)   ■  oT  y  sec(x)   csc(x)   esc    (y)  -  tan(x) 

(A-9) 

and, 
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12  2 

„  -z  esc   (x)    {y   esc   (y)-cot(y)}+cot(x)-cot(i) 

rl    ^  c • 

2d       y{l-2cot2(i)}+cot(i){cot(x)-cot(y)}+cot(x)cot(y) 


(A-10) 


When  a  force  F     is   inserted,    the  same   relations 
P 

can  be  used  to  define  a  new  expression  for  the  stabil- 
ity  number. 

Summation  of  forces   in  the  direction  of   the  x  axis 
(Figure   3-11)    gives 


F      cos    (CEO)  +  C     cos(x)   -  Psin(v)  (A-ll) 

P  r 


Summation   of   moments   around  point   0   gives: 


P{Rsin(<J>)}   +  Cr   (00')  +  F      (OG)   -  w     -   0 


(A-12) 


P   is   obtained   from  Equation    (A-ll): 


F      cos (CEO)   +  C     cos(x) 

=  -2 a    t  \     <A_13> 

sin(v) 


Substitution  of  Equations  (A-5),  (A-6),  and  (A-13)  in 
Equation  (A-12)  results  in: 

F  cos(CE0)+C  cos(x) 

{-2 ;  ,  >r }Rsin(<{.)+C  {Rsin(4.)csc(u-v)} 

sinv.v)  r 

+F  (OG)-^-  E  -  0  (A-14) 

P         Li. 
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The   new  stability  number  for  toe  failure  is 
obtained  by   substituting  Equations    (A-l)   and   (A-3)    into 
this   equation  and  simplifying: 

12   F 
E ^  {C°^)}  f  csc(x)csc(y)sin(  tf+OG} 


c  6csc    (x)csc(y)sin(  $){— - (-csc(u-v)cos(x-u)} 

(A-15) 

with   the  angles  u  and  v  obtained  from  Equations  A-8  and 
A-9. 

If   the  failure  surface  extends   below  the  toe   of 
the   slope: 

3  3 

Wd  =  -j|-  E  +  -2|-  {2n2  -  2n  sin(<|>)   csc(x)   csc(y)}    (A-16) 

where  r  defines  the  distance  from  the  toe  of  the  slope 
to  the  point  of  the  rupture  surface  on  the  ground  as 
AA'  =  rw  (Figure  3-10  in  Taylor,  1937): 


n  =  j   {cot(x)-cot(y)-cot(i)+sin(<{i)  csc(x)  csc(y)}(A-17) 


With  this  expression  of  Wd  the  following  stability 
number  is  obtained  for  failure  below  the  toe: 


,  12F 

{E+6n  -6nsin($)csc(x)csc(y)} -2-  A 

Ca  =  yer 

c     6csc  (x)csc(y)sin(  ♦)lsl„/ y\ 'csc(u-v)cos(x-u)  \ 


(A-18) 


143 


where 


A  .  cob(CEO)   H  C8c(x)c8c(y)8ln(  fl   +  OG  (A-19) 

sin(v)      2 


If   only   a  fraction   of   F      is   assumed  mobilized,    F 

'  p  P 

is  divided  by  the  proper  number  in  Equations  (A-15)  and 
(A-18). 
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APPENDIX  B 


Determination  of  the  Pile  Length  Above  the  Failure  Surface 

Equations  are  derived  to  express  the  pile  length 
above  the  critical  surface,  CE,  in  terms  of  the  angles 
x  and  y  defining  the  critical  circle  (Figure  9). 
Hence,  each  time  a  new  potential  critical  surface  is 
considered  (defined  by  x  and  y),  a  new  unit  force,  F  , 
function  of  CE,  is  determined  and  used  in  the  stability 
computations  for  that  particular  surface. 

Figures  34  and  35  show  the  geometry  of  the  slope. 
S'  is  the  horizontal  distance  between  the  toe  and  the 
center  of  the  circle,  and  S  is  the  horizontal  distance 
between  the  toe  and  the  pile.  CE  is  the  pile  length 
above  the  failure  surface.   OG  is  the  moment  arm  of  the 
force  F  .  x  is  the  angle  between  the  chord,  AB,  and 
the  horizontal,  and  2y  is  the  angle  AOB  of  the  circle. 
The  distance  S  is  an  input  parameter  since  it  defines 
the  location  of  the  pile  upslope. 
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Figure  34     Toe  Failure. 


l'»6 


Figure   35      Failure  below  the  Toe. 
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For  a  toe  failure,  the  geometry  of  Figure  34  is 
considered.  The  following  geometric  relationships  exist 
between  the  different  segments  and  angles  in  this  fig- 


ure: 


AC £p-  (B-D 

cos(i) 


OAB  =  OBA  =  j-y  (B-2) 


S'  =  R  cos  (OAB  +  x)  (B-3) 


OAC  =  OAB  -  (i-x)  (B-4) 


2    2     2 
-1  ,0C  -  R  -  AC  ,         (      Sv 
OCA  =  cos    { _2AC(0C) }         (B  5) 


and 


OCE  =  OCA  +  90-i  (B-6) 

Hence,  using  Equations  (B-l)  to  (B-6),  CE  is  obtained 
as: 


CE 


2(0C)cos(0CE)  +  \|{2(OC)cos(OCE)}2-4(OC2-R2)  (B_?) 


2 


If   the  pile  is  at  the  left  of   the  origin,   0,    the 
angle  OCE   is    replaced  in  this   equation   by: 


OCE  =  |  Ti  +  i  -  OCA  (B-8) 
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For  cases  where  the  failure  surface  passes  below 
the  toe,  the  geometry  of  Figure  35  Is  considered.   The 
following  geometric  relations  exist: 

AA'  =  n  H  (B-9) 

where  n  =  —   {cot(x)-cot(y)-cot(i)+sin(  <J>)  csc(x)  csc(y)} 

(B-10) 

AC  =  D/cos(i)  (B-ll) 

OAA'  =  OAB  +  x  (B-12) 


and 


-2   R2  -  OA'2  -  AA'2 

0A'A  ■ cos    i^^mhoPT  }       (B~13) 


OA'C  =  n  -  OA'A  -  i  (B-14) 


nrA'      ~2  f0A/2  -  A'C2  -  PC2,       ,_  ... 
OCA  =  cos    { _2  A>c  (QC) }      (B-15) 


OCE  =  OCA'  +  (y  -  i)  (B-16) 


CE  is  then  calculated  by  introducing  the  values 
obtained  by  Equations  (B-8)  to  (B-16)  in  Equation  (B- 
7).   If  the  pile  is  at  the  left  of  the  origin,  the 
angle  OCE  is  given  by  Equation  (B-8). 
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The  force  F  is  assumed  to  act  at  the  centroid  of 
P 

the  pressure  diagram  in  a  direction  parallel  to  the 
tangent  to  the  circle  at  the  point  of  intersection 
between  the  pile  and  the  failure  surface.  Hence,  the 
moment  arm,  OG,  of  this  force  is: 

OG  -  R  -  G'E  cos  (CEO)  (B-17) 

where  G'E  is  the  distance  from  the  centroid  of  the 
pressure  diagram  to  the  failure  surface  and  the  angle 
CEO  is  given  by: 


150 


APPENDIX  C 

Slope  Stability  Program 

After  reading  the  input  data,  the  slope  stability 
program  performs  the  following  steps. 

1)  Initializes  an  angle  of  friction  ij> 

2)  Computes  the  stability  number  for  every  possible 
surface  through  the  slope 

If  piles  are  used,  the  following  additional 
steps  are  necessary  for  the  calculation  of  a 
stability  number 

a)  The  distance,  measured  along  the  length  of 
the  pile  from  the  face  of  the  slope  to  the 
surface  under  investigation,  is  computed 

b)  The  reaction, from  the  piles  on  the  slope, 
is  calculated 

c)  The  moment  arm  from  the  center  of  the  phi- 
circle      is   found 

d)  The  stability  number  for  the  surface  is 
calculated 

3)  Finds  the  critical  stability  number  after  all 
possible  surfaces  have  been  investigated 

4)  Uses  the  critical  stability  number  to  calculate 
a  factor  of  safety  with  respect  to  cohesion 

5)  Checks  for  convergence 

If  the  difference  between  the  factor  of  safety 
with  respect  to  cohesion  and  the  one  with 
respect  to  friction  is  larger  than  the 
tolerance  the  program  repeats  steps  1  to  4. 
If  the  difference  is  less  than  the  tolerance, 
it  exits  the  itteration 

6)  Prints  out  the  input  and  output 
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The  following  data  is  required  for  the  design 
program.  All  input  is  free  formated. 


CARD  1 


slope 

h 

s 


CARD  2  pilein 
nosur 


CARD  3  phiav 

c 
gamma 

CARD  4  fphi 

sfdif 


DATA  CARDS 

slope  angle  in  degrees 

height  of  slope 

horizontal  distance  from  the  toe 

of  the  slope  to  the  piles 

indicates  the  presence  of  piles 
indicates  the  type  of  analysis 
needed.  For  pilein=0  and 
nosur=0,  the  unreinforce  slope 
is  analyzed.  For  pilein=l  and 
nosur=0,  the  reinforced  slope 
is  analyzed  and  the  new  critical 
surface  is  found. 
For  pilein=0  and  nosur=l,  the 
reinforced  slope  is  analyzed 
assuming  that  the  piles  did  not 
alter  the  original  critical  surface 

friction  angle  for  the  soil,  in 

degrees 

cohesion  intercept  of  the  soil 

unit  weight  of  the  soil 

initial  factor  of  safety  with  respect 
to  friction;  usually  taken  as  1.0 
the  difference  between  the  safety 
factor  with  respect  to  friction  and 
the  one  with  respect  to  cohesion, 
used  when  checking  for  convergence. 
A  difference  of  0.001  is  sufficient 


CARD  5  dO 
dl 


CARD  6  fpar 


the  pile  diameter, b. 

the  center  to  center  distance  between 

the  piles 

degree  cf  mobilization  of  Fp.  For 
fpar=l,  the  total  force,  Fp,  will 
be  used.  For  Fp=0  the  force  will  be 
divided  by  the  factor  of  safety 
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c 

c 

c DATA  INPUT 


c 
c 


c 

c   slope    =  slope  angle  in  degrees 

c   h        =  slope  height 

c   s       -   the  horizontal  distance  from  the  toe  of 

c  the  slope  to  the  piles 

c   pilein   =  The  two  quantities  are  used  together  to  indi- 

c   nosur      cate  the  presence  of  the  pile  and  the  type 

c  of  analysis  needed.  When  no  piles  are 

c  present,  both  pilein"0,  and  nosur=0 

c  When  piles  are  present  and  they  change  the 

c  critical  surface,  pilein»l,  and  nosur=0 

c  When  piles  are  present  but  are  not  expected 

c  to  change  the  critical  surface,  pilein=0,  and 

c  nosur=l 

c   phiav    ■  friction  angle  of  the  soil  in  degrees 

c   c       ■  cohesion  of  the  soil 

c   gamma    ■  unit  weight  of  the  soil 

c   fphi     =  initial  factor  of  safety  with  respect  to 

c  friction 

c   sdif     ■  difference  between  the  safety  factor  with 

c  respect  to  friction  and  the  one  with  respect 

c  to  cohesion  that  will  be  used  to  check 

c  convergence.     A  difference   of  0.001  will  be 

c  sufficient 

c       dO  =  pile   diameter 

c       dl  =  center  to  center  distance  between  the  piles 

c       fpar  =  0  if  Fp  is  partially  mobilized.     Otherwise, 

c  it   is   equal  to   1 

c       dfac  =  a  parameter  in  subroutine  phycir.     The  max 

c  depth  of   the  failure  surface  measured  from 

c  the  top  of   the  slope  is    limited  to   dfac*h 

c 

c ; 1 — 

c 

c 

c         MAIN  PROGRAM 

parameter (pi=3 . 14 159 ) 

common  /all/  slope, h,cosslp,cgch,oc, phi, xmax,ymax, 
1  nmax,dmax,tanslp, 

1  snmax, m, cots lp , sinphi , pilein, gamma , s , cosx, t anx , 

2  d,n,x,y,c8cx,cscy,ge,fp,ce,ceo,og,r,oach,direc 
common  /pres/  d0,dl,d2,tanphi,c,p(2), phiav 

c 

common   /phi/   cemax,fpmax,rmax,f c 
common   /pc/  ii,jj,kk, 11, nosur, fpar 


real  nmax 
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integer  pilein.fpar 

c 

read  (5,*), slope, h, s 
read  (5,*),pilein,nosur 
read  (5,*),phiav,c,gamma 
read  (5,*),d0,dl 
read  (5,*),fpar 

c 

d2=dl-d0 

c 

c  In  defining  the  different  failure  surfaces, 

c  x  goes  from  15  to  m  degrees  and  y  goes 

c  from  25  to  50  degrees 

c 

m^nintCslope)-! 
ii=25 
jj=50 
kk=15 
ll=m 
fc=1.0 

c 

phiav=phiav*pi / 180. 

slope=slope*pi/180. 

tanphi=tan(phiav) 

cosslp=cos (slope) 

cotslp=l./tan(slope) 

tanslp=tan(slope) 

c 

c  a  new  critical  surface  is  found  after  the  addition 

c  of  the  piles 
c 

c  initialize  a  friction  angle 
c 
c 
100   phi  =atan(tan(phiav)/fphi) 

c 

sinphi=sin(phi) 

c 

c  find  the  critical  stability  number 


c 
c 
c: 

c 


call  phicir 


if  (snmax.le.0.0)  then 

print*, 'You  are  overdesigning  the  slope' 

stop 

endif 


c 
c: 


f c=c / ( s nmax*gamma *h ) 
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c  when  phi  is  equal  to  zero 
c 

if  (phiav.le.0.0)  then 

f  3=f  C 

go  to  05 
endif 
c 

c  check  for  convergence  of  the  safety  factors 
c 

diff=abs(fc-fphi) 
if  (diff.gt.sfdif)  then 
fphi=(fphi+fc)/2. 
go  to  100 
else 

fs-(fphi+fc)/2. 
endif 
05     continue 
c 

c  one  well  defined  surface  that  does  not  change  with  the 
c  presense  of  the  piles 
c 

if  (nosur.eq.l.and.pilein.eq.O)  then 

c 

pilein=l 
c 

ii-ymax*180/pi 
jj=»ymax*180/pi 
kk=«xmax*180/pi 
ll=xmax*180/pi 
c 

fphi=l. 
fc-1. 
c 

c  intialize  a  friction  angle 
c 
101    phi=atan(tan(phiav)/fphi) 
sinphi=sin(phi) 
c 

c  find  the  critical  stability  number 
c 

call  phicir 

c: : : : : 

c 

if  (snmax.le.0.0)  then 

print*, 'You  are  overdesigning  the  slope' 

stop 

endif 
c 
c: : : : : 

f  c°c / ( s  nma x*gamma  *h ) 
c 
c  when  phi  is  equal  to  zero 
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c 

if  (phiav.le.O.O)  then 

fs-fc 

go  to  06 

endif 
c 
c  check  for  convergence  of  the  safety  factors 

c 

diff=abs(fc-fphi) 

if  (diff.gt.sfdif)  then 
fphi=(fphi+fc)/2. 
go  to  101 
else 
fs=(fphi+fc)/2. 
endif 
06      continue 
endif 
c 

slope=slope*180. /pi 
phiav=phiav*180. /pi 
xmax=xmax*180./pi 
ymax=ymax* 1 80 . /p  i 
c 
c 

write  (6,440)  slope, h,phiav,c, gamma 

440  format  ('SLOPE  CONFIGURATION'// 
1  5x, 'Slope  Angle  -  ',f5.1/ 

1  5x, 'Slope  Height  -  ',f5.1// 

1  'SOIL  PARAMETERS'// 

1  5x, 'Friction  Angle  -  ',f5.1/ 

1  5x, 'Cohesion  -  ',f6.1/ 

1  5x,'Unit  Weight  =  ',f6.1//) 

c 

if  (pilein.eq.O.and.nosur.eq.0)  go  to  450 

write  (6,441)  dO.dl.s 

441  format  ('PILE  DIMENSIONS'// 
1  5x, 'Pile  Diameter  -  ',f4.1/ 

1  5x, 'Center-to-Center  Distance  ■  ',f5.1/ 

1  5x, 'Distance  Upslope  =  ',f5.1//) 

c 

if    (pilein.eq.l.and.nosur.eq.O)   then 

write    (6,442) 

442  format    ('GEOMETRY  OF  THE  CRITICAL   SURFACE'// 

1   5x,  "The   critical  surface  of   the  slope   changes'/ 
1   5x,'with  the  insertion  of   the  piles'//) 

else  if    (pilein.eq.O.and.nosur.eq.l)   then 
write   (6,443) 

443  format    ('GEOMETRY  OF  THE   CRITICAL   SURFACE'// 

1   5x,'The   critical  surface  of   the  slope   does  not   change'/ 

1   5x,'by  the  presence   of   the  piles'//) 

endif 

c 

write   (6,444)  xmax,ymax 
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444  format  (5x,'The  angles  that  define  the  circle  are:'/ 
1  lOx, 'x-  ',f5.1/ 

1  10x,'y«  ',f5.1//) 
c 

if  (nmax.lt.O. )  then 
write  (6,445) 

445  format    (5x, 'Toe  failure'//) 

else  if    (nmax.gt.O)   then 
hnmax=nmax*h 
write   (6,446)  hnmax 

446  format    (5x, 'Failure  below  the  toe.       The  critical  surface'/ 
1   5x, 'intercepts    the   ground  surface    ',f7.2, 'units    left    '/ 

1   5x,'of   the  toe'//) 
endif 
c 

if    (fpar.eq.l)   then 
write    (6,447) 

447  format    ('UNIT  REACTION   FORCE  ON  THE  SLOPE  EXERTED   BY  THE  PILES'// 
1   5x, 'Fp   is   assumed  to  be  totally  mobilized'//) 

else  if    (fpar.ne.l)   then 
write   (6,448) 

448  f ormat  '  ( 'UNIT  REACTION   FORCE  ON  THE   SLOPE  EXERTED   BY  THE  PILES'// 
1   5x, 'Fp   is   assumed  to  be  partially  mobilized,'/ 

1   5x, 'and  is   equal  to  Fp/FS'//) 
endif 
c 
c 

write    (6,449)   fpmax,cemax 

449  format    (5x, 'The  unit   reaction  force  on  the  slope  provided'/ 
1   5x,'by   the  piles   is    ',f7.0// 

1   5x, 'The  distance  from  the  ground  to  the  critical  surface,'/ 
1   5x,'CE,    is    'f5.1//) 
c 

450  write  (6,451)  sfdif,  fs 

451  format  ('FACTOR  OF  SAFETY'// 

1  5x, 'The  tolerable  difference  between  Fc  and  Fphi'/ 

1  5x, 'is  ',f7.4// 

1  5x, 'The  safety  factor  of  the  slope  is   ',f5.2) 


c 
c 


c 
c 

c- 


8  top 

end 


subroutine  phicir 

c 

c 

c  subroutine  that  computes  the  critical  surface  and 

c  stability  number 

c 

real  n,nmax,nfact 

common   /all/   s lope, h,cosslp,cgch,oc, phi, xmax,ymax, 
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c 
c 
c 


c 
c 


c 
c 


c 
c 


1   nmax,dmax,tanslp, 

1  snmax,m,cotslp,sinphi,pilein, gamma, s,cosx, tanx, 

2  d,n,x,y,cscx,cscy,ge,fp,ce,ceo,og,r,oach,direc 

common    /phi/   cemax,fpmax, rmax,f c 
common    /pc/   ii, j j,kk,ll,nosur,fpar 

integer  m, pilein, fpar 

parameter (pi=3. 14 159, dfac=8. 50, nfact=169.) 


if    (fc.le.1.0)    f-1.0 
if    (fc.gt.1.0)   f=fc 
if    (fpar.eq.l)   f-1.0 


snmax=0.0 

do  20  numy=ii,jj 
do  10  numx=kk,ll 


x=numx*pi/180. 
y=numy*pi/180. 
cosx=cos(x) 
tanx=tan(x) 
cotx=l . /tan(x) 
coty=l./tan(y) 
cscx=l . /sin(x) 
cscy=l ./sin(y) 
secx=l . /cos(x) 
cosx=cos(x) 
siny=sin(y ) 

d=0.5*(cscx*cscy-cotx*coty+l ) 
n=0.5*(cotx-coty-cotslp+sinphi*cscx*cscy) 


if  (n.le.0.0)  then 
c. 

c  the  failure  surface  passes   through   the  toe 
c 

if  (y.gt.x.and.d.ge.dfac)  go  to  20 


c 
c 


if  (phi.le.O.O)go  to  06 

h2dn=0.5*cscx**2*(y*cscy**2-coty)+cotx-cotslp 
h2dd=l./3.*(l.-2.*cotslp**2)+cotslp*(cotx-coty)+cotx*coty 
h2d  =h2dn/h2dd 

cotu=h2d*y*secx*cscx*cscy**2-tan(x) 
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u=atan(l./cotu) 
8lnuv=h2d*sin(u)*cscx*C8cy*sinphi 
C8cuv=l . /sinuv 

angl=asin(sinuv) 

v=u-angl 

sinv=8in(v) 


c 


06    if  (pilein.le.O)  then 
c 

c  ***  without  piles   *** 
c 

if  (phi. gt. 0.0)  then 
cfwhn=h2dn 
cf whd=2 . *cotx/tan( v)+2 . 

else  if    (phi. le. 0.0)   then 
cfwhn=(l.-2*cotslp**2)/3.+cotx*coty+cotslp*(cotx-coty) 
cfwhd=2.*y*c8cx**2*cscy**2 
endif 
c 

sn=cfwhn/cfwhd 
else  if    (pilein.ne.O)   then 
c 

c     ***     with  piles     *** 
c 

call   momarm 
c  =========== 

eq4=l.-(2.*cotslp**2)+(3.*cotslp*cotx)-(3.*cotslp*coty) 
1   +(3.*cotx*coty) 
c 

if    (phi. gt. 0.0)   then 
cfwhn=eq4-((12.*fp)/(f *gamma*h**3)*(cos(ceo)/sinv*h/2.* 
1   cscx*cscy*sinphi+og)) 

cfwhd=6.*cscx**2*cscy**2*siny*sinphi*(cosx/sinv+ 
1   cscuv) 

else  if  (phi. le. 0.0)  then 
cfwhn=eq4-12.*fp*og/(f *gamma*h**3) 
cfwhd=6.*y*c8cx**2*cscy**2 
endif 
c 

sn=cfwhn/cfwhd 
endif 
c 

else 

c 

c  the  failure  surface  passes  below  the  toe 

c 


c 
c 


if  (n.gt.nfact)  go  to  20 

if  (y.gt.x.and.d.ge.df ac)  go  to  20 


if  (phi. le. 0.0)  go  to  07 
h2dn=0 . 5*cscx**2* (y *cs  cy**2-coty )+cotx-cotslp-2 . *n 


159 


h2dd-l./3.*(l.-2.*cotslp**2)+cotslp*(cotx-coty)+cotx*coty 
1  +2.*n**2-2.*n*sinphi*cscx*cscy 
h2d=h2dn/h2dd 

cotu=h2d*y*secx*cscx*cscy**2-tan(x) 

u=atan(l./cotu) 
sinuv=h2d*sin(u)*cscx*cscy*sinphi 

cs  cu v= 1 . /s  i  nuv 
angl=asin(sinuv) 

v=u-angl 

sinv=sin(v) 

07  if    (pilein.eq.O)    then 


c 

c  ***  without  piles   *** 

c 

if  (phi. gt. 0.0)  then 

cfwhn=h2dn 
cfwhd=2.*cotx/tan(v)+2. 

else  if  (phi. le. 0.0)  then 
cfwhn=(l.-2.*cotslp**2)/3.+cotx*coty+cotslp*(cotx-coty)+2*n**2 

cfwhd=2.*y*cscx**2*cscy**2 
endif 

sn=cfwhn/cfwhd 

else  if  (pilein.ne.O)  then 
c 

c  ***  with  piles   *** 
c 

call  momarm 


c 
c 

c: 


if  (oach.gt.r)  go  to  20 
eq4=l.-(2.*cotslp**2)+(3.*cotslp*cotx)-(3.*cotslp*coty) 

1  +(3.*cotx*coty) 

if  (phi. gt. 0.0)  then 

cf whn=(eq4+3 . *(2 . *n**2-2 . *n*sinphi*cscx*cscy ) )-12 . *f p / 
1  (f*gamma*h**3)*(cos(ceo)/sinv*h/2.*cscx*cscy*sinphi+og) 

cfwhd=6.*cscx**2*cscy**2*siny*sinphi*(cosx/sinv+cscuv) 

else  if  (phi. le. 0.0)  then 
cfwhn=eq4+6*n**2-12.*fp*og/(f*gamma*h**3) 

cfwhd=6.*y*cscx**2*cscy**2 

endif 
sn=cfwhn/cfwhd 
endif 

endif 


if    (sn.gt.snmax)    then 

snmax=sn 

ymax=y 

xmax=x 
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c 

c 


dmax=d 

nmax=n 

cemax=ce 

fpmax*fp 

rmax=r 


p  •  •  •  • 

•  ••••*■••• 

c 

c 

c 

endif 

10 

continue 

20 

continue 

c 

return 

end 

subroutine   momarm 


c 

c 

c  subroutine  that  calculates  the  moment  arm  of 

c  force  Fp 

c 

common  /all/  slope, h,cosslp,cgch,oc, phi, xmax,ymax, 

1  nmax,dmax,tanslp, 

1  snmax.m, cots lp,sinphi,pilein, gamma, s.cosx.tanx, 

2  d,n,x,y,cscx,cscy,ge,fp,ce,ceo,og,r,oach,direc 
real  n 

parameter  (pi=3. 14159) 


r=h/2.*cscx*cscy 
oab=pi/2.-y 

if  (n. It. 0.0)  then 

the  failure  surface  passes  through  the  toe 


c. 
c 

c, 


ac=s/cosslp 

oabx=oab+x 

sch=r*cos(oabx) 

oa  c=oa  b-s lop  e+x 

oc=sqrt(r**2+ac**2-2.*r*ac*cos(oac)) 

oca=acos((r**2-oc**2-ac**2)/((-2.)*ac*oc)) 


if    (sch.gt.s)  then 
pile  to  the   left   of   origin 
oce=3 . /2 . *pi+s lope-oca 
else  if    (sch.lt.s)   then 
pile   to  the   right   of   origin 
oce=pi/2.-slope+oca 
else  if    (sch.eq.s)   then 
pile  underneath   the   origin 
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c 

c 


ce=*r-oc 
og=oc+cgch 
go  to  111 
endif 

else   if    (n.gt.0.0)   then 


c 

c        the     failure     surface     passes      below     the     toe 

c 

aach=n*h 

achc=s/cosslp 

oaach=oab+x 

oach=sqrt  (aach**2+r  **2-2 .  *r  *aach*cos  (oaach ) ) 

oacha=acos((r**2-oach**2-aach**2)/((-2.)*aach*oach)) 

oachc=pi-slope-oacha 

oc^sqrt (oach**2+achc**2-2 . *oach*achc*cos (oachc ) ) 

ocach=acos((oach**2-achc**2-oc**2)/((-2.)*achc*oc)) 

sch=oach*cos(pi-oacha) 

c 

if    (sch.gt.s)   then 
c       pile   on   left   of   origin 

oce=3 ./2. *pi+s lope-ocach 

else   if(sch.lt.s)    then 
c       pile  to  the   right   of   origin 

oce=ocach+pi/2.-slope 

else   if    (sch.eq.s)   then 
c       pile  underneath   origin 

ce=r-oc 

og=oc+cgch 

go   to   111 

endif 

endif 


ce=(2.*oc*cos(oce)+sqrt((2*oc*cos(oce))**2-4.*(oc**2 
1   -r**2)))/2. 
cosc=cos(ceo) 

cosc=(oc**2-ce**2-r**2)/((-2.)*ce*r) 

ceo=acos(cosc) 

call  presur 

gche=ce-cgch 
eg=gche*cosc 
if    (direc.eq.l)   then 
og=r-eg 
111      continue 


c 
c 


return 
end 


c 
c 
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c 

subroutine  presur 

c 

c  subroutine  that  calcuates  pressure  Fp 

c 

parameter  (pi=3. 14159) 

common   /all/   s lope, h,cosslp,cgch,oc, phi, xmax,ymax, 
1   nmax,dmax,tanslp, 

1  snmax.m, cotslp, sinphi, pilein, gamma, s, cos x, tanx, 

2  d,n,x,y,cscx,cscy,ge,fp,ce,ceo,og,r,oach,direc 
c 

common   /pres/  d0,dl,d2,tanphi,c,p(2) ,phiav 

real  nphi,nphil,nphi2 
c 
c 

c  equation  in  the  form  p=a+bz   gives    the  pressure   on  the  pile 

c  as   a   function   of   depth  z 

c 

z-0.0 

do   10   1=1,2 

nphi=(tan(pi/4.+phiav/2.))**2 

nphil=(nphi)**(.5) 

nphi2=(nphi)**(-.5) 

parl=(dl-d2)/d2*nphi*tan(pi/8.+phiav/4.)*tanphi 

par2=(2.*tanphi+2*nphil+nphi2)/(nphil*tanphi+nphi-l) 

par3=nphil*tanphi+nphi-l 
c 

if  (c.le.0.0)  then 
c 

c      cohesionless  soil 
c 

p(i)=gamma*z/nphi*(dl*(dl/d2)**par3*exp((dl-d2)/d2*nphi* 
1  tanphi*tan(pi/84phiav/4))-d2) 

else  if  (phiav.le.0.0)  then 
c 

c  cohesive     soil 

c 

p(i)°c*(dl*(3*alogl0(dl/d2)+(dl-d2)/d2*tan(pi/8)) 
1  -2*(dl-d2))+gamma*z*(dl-d2) 

else 
c 

c  c-phi     soil 

c 

p(i)=c*dl*(dl/d2)**par3*(l/(nphi*tanphi)*(exp(parl)-2*nphil 

1  *tanphi-l)+par2)-c*(dl*par2-2*d2*nphi2)+gamma*z/nphi* 

2  (dl*(dl/d2)**par3*exp(parl)-d2) 
endif 

z»ce 
10  continue 

c 

c  to  find  the  force    (ie.the  area  of   the  pressure  diagram 

c  and  its    center  of   action) 
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c 
c 


al"ce*p(l) 
a2=ce*(p(2)-p(l))/2. 

yl»ce/2. 

y2=2*ce/3. 

alyl=al*yl 

a2y2=a2*y2 

sum=alyl+a2y2 

ftot=al+a2 

y3=sum/f tot 

cgch=y3 

fp=ftot/dl 


return 
end 
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DATA 


30.      45.        25. 

1   0 

10.      500.      125. 

1.      .001 

3.        7.5 

0 
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RESULTS 

SLOPE  CONFIGURATION 

Slope  Angle  =  30.0 
Slope  Height  =  45.0 

SOIL  PARAMETERS 

Friction  Angle  =  10.0 
Cohesion  =  500.0 
Unit  Weight  =  125.0 

PILE  DIMENSIONS 

Pile  Diameter  =  3.0 
Center-to-Center  Distance  =   7.5 
Distance  Upslope  =  25.0 

GEOMETRY  OF  THE  CRITICAL  SURFACE 

The  critical  surface  of  the  slope  changes 
with  the  insertion  of  the  piles 


The  angles  that  define  the  circle  are: 
x=  23.0 
y=  35.0 


Toe  failure 


UNIT  REACTION  FORCE  ON  THE  SLOPE  EXERTED  BY  THE  PILES 

Fp  is  assumed  to  be  partially  mobilized, 
and  is  equal  to  Fp/FS 


The  unit  reaction  force  on  the  slope  provided 
by  the  piles  is  21424. 

The  distance  from  the  ground  to  the  critical  surface, 
CE,  is   16.5 
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FACTOR  OF  SAFETY 

The  tolerable  difference  between  Fc  and  Fphl 
is  0.0010 

The  safety  factor  of  the  slope  is   1.33 
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APPENDIX  D 

Pile  Analysis  Program 

After  reading  the  input  data,  the  pile  analysis 
program  performs  the  following  steps. 

1)  Computes  the  force  F  on  the  pile 

2)  Constructs  an  array  and  solves  the  simultaneous 
equations  produced  by  the  finite  difference 
method.  The  displacement  of  the  pile  below 
the  critical  surface  is  thus  obtained 

3)  Calculates  the  slope,  bending  moment  and  shear 
force  on  the  pile  section  below  the  critical 
surface 

4)  Calculates  the  displacement,  slope  bending 
moment,  and  shear  force  for  the  pile  section 
above  the  failure  surface,  using  a  closed 
form  solution 

5)  Prints  out  the  input  and  output 
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The  following  data  Is  required  for  the  pile  analysis 
program.  All  input  is  free  formated. 


CARD  1   phiav 


c 

gamma 


CARD  2  CE 
EL 
EB 
BD 


CARD 

3 

MT 

CARD 

4 

dO 
dl 

EI 

CARD 

5 

ksO 

ksl 


krO 


krl 


CARD  6   nl 


DATA  CARDS 

friction  angle  for  the  soil,  in 

degrees 

cohesion  intercept  for  the  soil 

unit  weight  of  the  soil 

depth  from  the  slope  face  to  the 

critical  surface 

depth  of  pile  below  the  critical 

surface 

depth  of  pile  from  the  critical 

surface  to  the  bedrock 

depth  of  pile  embedded  in  the 

bedrock.  If  bedrock  is  not 

present,  the  quantities  EB  and  BD 

can  be  used  to  indicate  the  depth 

of  two  different  layers  of  material 

below  the  critical  surface.   If 

the  material  below  the  critical 

surface  is  uniform,  one  of  the 

two  quantities  can  be  set  equal 

to  zero. 

number  of   discretized  points   below 
the   critical  surface 

the  pile   diameter,    b. 

the   center-to-center   distance   between 

the  piles 

pile  stiffness 


the  coeff 
subgrade 
surface 
the  coeff 
subgrade 
layer  EB 
the  coeff 
subgrade 
of  layer 
the  coeff 
subgrade 
the  pile 


icient  of  the  horizontal 
reaction  at  the  critical 

icient  of  the  horizontal 
reaction  at  the  end  of 

icient  of  the  horizontal 

reaction  at  the  beginning 

BD 

icient  of  the  horizontal 

reaction  at  the  end  of 


index  greater  or  equal  to  zero. 
Indicates  the  variation  of  ks  with 
depth,  for  the  layer  EB.   For 
example,  if  nl=0,  ks  is  constant 
with  depth.   If  nl*!,  ks  increases 
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linearly  with  depth. 
n2       same  Index  as  nl.   Indicates  the 
variation  of  ks  with  depth  for 
layer  BD 

CARD  7  be       boundary  condition  at  the  top  of 
the  pile,  equals  to  1,  2,  3,  or  4 
(see  Chapter  4) 

hr      Indicates  the  degree  of  hardness 
of  layer  BE.  For  hr»l,  the  layer 
is  composed  of  hard  rock  (stiffer 
than  the  pile  itself).  For  hr=0, 
this  layer  is  composed  of  soft  rock 
or  soil 

kopt     gives  the  user  the  option  to  input 
ks  for  every  single  point  of  the 
discretized  pile  below  the  critical 
surface.  For  kopt=0,  a  total 
number  of  MT  additional  data  cards 
(following  card  7)  are  required. 
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c 

c 

C DATA  INPUT     

c 

c ■ 

c 

c    phiav   =  friction  angle  of  the  soil  in  degrees 

c    c       =  cohesion  of  the  soil 

c    gamma   =  unit  weight  of  the  soil 

c    CE      ■  depth  from  the  slope  face  to  the  critical 

c  surface 

c    EL      =  depth  of  the  pile  below  the  critical  surface 

c    EB      =  distance  between  the  critical  surface  and 

c  the  bedrock 

c  BD  =   depth    of    the    pile    below    the    bedrock 

c  mt  =   number   of   points    the  pile  will   be    descretized 

c  to 

c  dO  =   pile    diameter 

c  dl  =    center-to-center   distance    between   the   piles 

c  ksO  =   coefficient    of    subgrade    reaction  at    the 

c  critical   surface 

c  ksl  =   coefficient    of    subgrade    reaction   at    the    end 

c  of    layer   EB 

c    krO     =  coefficiet  of  subgrade  reaction  at  the  beginning 

c  of  layer  BD 

c    krl     =  coefficient  of  subgrade  reaction  at  the  end  of 

c  the  pile 

c    nl      =  index  that  indicates  the  variation  of  ks  with 

c  depth.   For  nl=0,  ks  is  constant,  for  nl-1 

c  ks  increases  linearly   (in  layer  EB) 

c    n2      =  same  as  nl,  for  layer  BD 

c    ei      =  pile  stiffneso 

c  be  =  boundary   condition   for   the  pile   top    (1,    2,    3,    or   4) 

c  hr  =  for  hr=l,    layer   BE    is    assumed   to   be   hard   rock 

c  for  hr=0,    it    is    assumed   to   be   soft    rock   or   soil 

c  kopt  =   if   kopt=0,    a   total  number   of   mt   additional   data 

c  carda    are    required   to  input    the    values    of    the 

c  subgrade    reaction  along   the   pile 

c 

c 

c 

c 

c  MAIN  PROGRAM 

dimension  a (85, 86) ,g(85,86) ,h(85,86) ,y(86) 

dimension  zl(86),z2(86),sl(86),s2(86) 

dimension  ml (86) ,m2(86) ,vl(86) ,v2(86) ,yl(86) 

dimension  ki(86) 
c 

parameter  (mtl=6,pi=3. 14159) 
c 

common  /all/  ce,d0,dl , phiav, c, gamma, p(2) 
c 

integer  be,  r,  hr 
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real  nl,n2,   ksO,ksl,ks,   krO.krl.kr,  ml,m2,   ki 

read(5,*)  phiav.c, gamma 
read(5,*)  ce,el,eb,bd 
read(5,*)  mt 
read(5,*)  dO.dl.el 
read(5,*)  ksO,ksl,krO,krl 
read(5,*)  nl,n2 
read(5,*)  bc,hr,kopt 

if  (kopt.le.O)  then 
do  480  i=3,mt+3 
read(5,*)  ki(i) 
480   continue 

endif 


c 
c 


write  (6,460)  phiav, c, gamma 

460  format  ('SOIL  PARAMETERS'// 

1    5x, 'Above   The   Critical    Surface'// 

1    lOx, 'Friction  Angle   =    ',f4.1,'    degrees'/ 

1    lOx, 'Cohesion   =   '  ,f 6.1/ 

1    lOx.'Unit  Weight   =    ',f6.1// 

1    5x,    'Below  The  Critical   Surface'//) 

if    (kopt.ne.O)    then 
write    (6,461)  ksO,ksl,krO,krl,kopt 

461  fornat    (lOx, 'Coefficients    of    Subgrade   Reaction'// 
1    20x,'ks0   =    ',e!0.3/ 

1  20x, 'ksl   =    ',el0.3/ 

1  20x,'kr0   =    ',el0.3/ 

1  20x,'krl   =   ',el0.3// 

1  10x,'kopt   =    ',12//) 

c 

else  if  (kopt.eq.O)  then 
write  (6,462) 

462  format  (lOx, 'Coefficients  of  Subgrade  Reaction'//) 

do  463  i=3,mt+3 
write  (6,464)  ki(i) 

464  format  (20x,el0.3) 

463  continue 

write  (6,465)  kopt 

465  format  (10x,/ 

1  10x,'kopt  =  ',12/ /) 

c 

endif 

c 

write  (6,466)  nl,n2 

466  format  (lOx, 'Variation  of  the  subgrade  reaction  coefficients  // 
1  20x,'nl  =  ',f3.1/ 

1  20x,'n2  =  ',f3.1//) 

c 

write  (6,467)  ce,el,eb,bd,dO,dl,ei,mt 

467  format  ('PILE  DIMENSIONS'// 
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1  5x, 'Length  of  pile  above  the  critical  surface  ■  ',f4.1/ 

1  5x, 'Length  below  the  critical  surface  =  ',f4.1/ 

1  lOx, 'Distance  between  the  critical  surface  and  the'/ 

1  lOx, 'bedrock  (or  hard  soil  layer)  =  ',f4.1/ 

1  lOx, 'Length  of  pile  embedded  in  the  bedrock  =  ',f4.1// 

1  5x, 'Pile  Diameter  =■  ',f4.1/ 

1  5x, 'Center-to-Center  Distance   =    ',f4.1/ 

1  5x, 'Stiffness    =   'el0.3// 

1  5x, 'Number   of   points    that   the  pile   length   below   the'/ 

1  5x, 'critical  surface  is  discretized  to  ■  ',i3//) 

if  (bc.eq.l)  print*, 'Boundary  Condition — Free  Head' 

if  (bc.eq.2)  print*, 'Boundary  Condition — Unrotated  Head' 

if  (bc.eq.3)  print*, 'Boundary  Condition — Hinged  Head' 

if  (bc.eq.4)  print*, 'Boundary  Condition — Fixed  Head' 


c 
c 
c 


phiav=phiav*pi/180. 
call  presur 
pl-p(l) 
P2=p(2) 
write  (6,469)  p(l),p(2) 

469  format  (5x,/ 

1   5x, 'Force  Diagram'/ 

1    lOx,  'force  on  the  pile  head  =  'f7.0/ 
1   lOx, 'force  acting  on  the  pile  at  the  critical'/ 
1    10X, 'surface  =  'f7.0//) 
c 

write  (6,470) 

470  format  (5x, 'Z' ,8x, 'DISPL(Y) ' , 7x, 'SLOPE(S) ',7x, 'MOMENT(M) ' , 
1    6x,'SHEAR(V)'//) 

c 

b=d0 
c 

ks=ksl-ks0 
kr=krl-kr0 
q=el/mt 
a3=q**4/ei 
fl=p2 

f2=(P2-pl)/ce 
c 

c  Finite  difference  solution  for  the  pile  section  below 
c  the  failure  surface  — solving  the  simultaneous  equations 

c*** ***************************************** ************* 

c 

c 

c  Constructing  the  array 

c 

c 

c     zero  array 

do  2  i=l,85 
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c 
c 
c 
c 


do  3  j-1,86 
a(i,j)=0. 
3  continue 

2  continue 

read  first  two  boundary  conditions 

if  (bc.eq.l.or.bc.eq.3)  then 

a(l,D=ce 
a(l,2)=2.*q-2.*ce 
a(l,3)=-4.*q 
a(l,4)=2.*q+2.*ce 

a((1i;mt+6)=2*q**3*(f2/(6*ei)*ce**3-fl/(2*ei)*ce**2) 

endif 

if  (bc.eq.l.or.bc.eq.2)  then 

a(2,l)=-l. 

a(2,2)=2. 

a(2,4)=-2. 

a(2;mt+6)=2*q**3*(fl/ei*ce-f2/(2*ei)*ce**2) 

endif 

if  (bc.eq.2.or.bc.eq.4)  then 

a(l,l)=-(ce**2) 
a(l,2)=2*ce**2-4*q*ce-2*q**2 

a(r,3)»S*q*ce 

a(l  4)=-2*ce**2-4*q*ce+2*q**2 

a(l^t+6)*4*q**3*(fl/(6*ei)*ce**3-f2/(24*ei)*ce**4) 

endif 

if  (bc.eq.3.or.bc.eq.4)  then 

a(2,l)=ce**3 

a(2  2)=6*q**2*ce+6*q*ce**2-2*ce**3 

a(2',3)  =  12.*q*(q**2-ce**2) 

a(2,4)=6*q*ce**2-6*q**2*ce+2*ce**3 

f(^t+6)=12*q**3*(f2/(120*ei)*ce**5-fl/(24*ei)*ce**4) 

endif 


c 
c 
c 


inputing  rest  of  array 

do  30  i=3,mt+3 

j=i-2 
a(i, j)=l- 

j=i-l 
a(i,j)="4. 

j=i 

m=i-3 
z=m*q 

if  (kopt.ge.l)  then 
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c 

c 


if  (z.le.eb)  then 

a(i, j)=6.+a3*b*(ks0+ks*(z/eb)**nl) 
elseif  (z.gt.eb)  then 

a(i, j)=6+a3*b*(krO+kr*((z-eb)/bd)**n2) 

endif 

elseif  (kopt.le.O)  then 
a(i, j)=6.+a3*b*ki(l) 


endif 

a(i,j)=-4. 

j=i+2 
a(i,j)=l. 
30  continue 

c 

c       inputing  last  boundary  conditions 
c 
c 

if  (hr.le.O)  then 
c  Soft  Rock,  Shear  and  Moment  equal  zero  at  EL 
c 
c 

i=mt+5 
a(i,mt+l)=-l. 
a(i,mt+2)=2. 
a(i,mt+4)=-2. 
a(i,mt+5)=l. 
c 

i=int+4 
a(i,iat+2)^l. 
a(i,mt+3)=-2. 
a(i,mt+4)=l. 
c 

elseif  (hr.ge.l)  then 
c  Hard  rock,  displacement  and  rotation  equal  zero  at  EL 
c 
c 

i=mt+5 
a(i,mt+l)=l. 
a(i,rat+2)=-4. 
a ( i, mt+3 )=a( mt+3, mt+3) 
a(i,mt+4)=-4. 
a(i,mt+5)=l. 
c 

i=mt+4 
a(i,mt+2)=-l. 
a(i,mt+4)=l. 
c 

i=mt+3 
a(i,mt+l)=0. 
a(i,mt+2)=0. 
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a(i,mt+3)=l. 
a(i,mt+4)=0. 
a(i,rct+5)=0. 
c 

endif 
do  70  i=l,mt+5 
do  80  j=l,mt+6 
c         print*,  a(i,j) 
80     continue 
70     continue 
c 

c 

c       calculating  the  auxiliary  quantities 

c 

c 

suml=0. 
sum2=0. 
do  111  l=l,mt+5 
do  222  k=l,mt+5 
if  (l.le.l)  go  to  12 
do  22  r-1,1-1 

suml=suml+g(k, l-r)*h(l-r, 1) 


22 

continue 

12 

g(k,l)=a(k,l)-suml 

suml=0. 

222 

continue 

c 

1-1 

do  333  j=i+l,at+6 

m=i 

if  (l.le.l)  go  to  13 

do  32  r=l,i-l 

sum2=sum2+g(i,i-r)*h(i-r, j) 

32 

continue 

13 

h(i, j)=(a(i, j)-sum2)/g(i,m) 

sum2=0. 

333 

continue 

111 

continue 

do  989  i=l,mt+5 
do  898  j=l,mt+6 
898   continue 
989   continue 
c 

c 

c      solving  for  the  unknown  displacements  y 

c 

c 

sum3=0. 

do  555  l=l,mt+5 

m=mt+6-l 
if  (l.le.l)  go  to  55 
do  65  r-1,-1-1 
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sum3=sum3+h(m,m+r )*y (m+r) 
65   continue 
55    y(m)=h(m,mt+6)-sum3 
sum3=0 . 

z»(m-3)*q 
if  (z.gt.el)  go  to  555 
555   continue 
c 

c 

c       calculating  slope  s,  moment  m,  and  shear  v  for 
c       z>0. 

c 

c 

do  666  m=3,mt+3 

z2(m)=(m-3)*q 
s2(m)=(y(m+l)-y(m-l))/(2*q) 
m2(m)=-ei*(y(m+l)-2*y(m)+y(m-l))/q**2 
v2(m)=-ei*(y(m+2)-2*y(m+l)+2*y(m-l)-y(m-2))/(2*q**3) 
666   continue 
c 
c ************************************************* ******** 

c  Close  form  solution  of  the  pile  section  above  the 

c  failure  surface 

c ********************************** *********************** 

c 

c 

c      calculating  deflection  y,  slope  s,  moment  m,  and 

c      shear  v  for  z<0 

c 

c 

a0=y(3) 

al=(y(4)-y(2))/(2*q) 
a2=(y(4)-2*y(3)+y(2))/(2*q**2) 
a3=(y(5)-2*y(4)+2*y(2)-y(l))/(12*q**3) 
c 

do  777  1=1, mt 1+1 
ql=ce/mtl 
zl(i)=-(i-l)*ql 
c 

yl(i)=a0+al*zl(i)+a2*zl(i)**2+a3*zl(i)**3+fl/(24*ei)*zl(i)**4+f2/ 
1  (120*ei)*zl(i)**5 

sl(i)=al+2*a2*zl(i)+3*a3*zl(i)**2+fl/(6*ei)*zl(i)**3+f2/(24*ei)* 
1  zl(i)**4 

ml(i)=-ei*(2*a2+6*a3*zl(i)+fl/(2*ei)*zl(i)**2+f2/(6*ei)*zl(i)**3) 
vl(i)=-ei*(6*a3+fl/ei*zl(i)+f2/(2*ei)*zl(i)**2) 
c 

777   continue 
c 

i=mtl+l 
442   write  (6,440)  zl(i),yl(i),sl(i),ml(i),vl(i) 
440   format  (2x,f 5. 1 ,5x,el0.3,5x,el0.3,5x,el0.3,5x,el0.3) 
i»i-l 

if  (i.lt.2)  go  to  441 
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c 


go  to  442 
441   continue 

do  443  m=3,mt+3 
write  (6,444)  z2(m) ,y(m) ,s2(m) ,m2(m),v2(m) 
444   format  (2x,f 5. 1 ,5x,el0.3 ,5x,el0.3,5x,el0.3,5x,el0.3) 

443   continue 

stop 
end 


c 


c 

c 
c- 

c- 


subroutine  presur 


c 

c  subroutine  that  calcultes  the  force  Fp 


c 


parameter  (pi=3. 14159) 

common  /all/  ce,dO,dl ,phiav,c, gamma, p(2) 


real  nphi.nphil ,nphi2 
c 
c 
c      equation  in  the  form  p=a+bz  gives  the  pressure  on  the  pile 

c      as  a  function  of  depth  z 

c 

d2=dl-d0 
tanphi=tan(phiav) 


c 

c 


z=0.0 

do  10  1-1,2 

nphi=(tan(pi/4.+phiav/2.))**2 

nphil=(nphi)**(.5) 

nphi2=(nphi)**(-.5) 
parl=(dl-d2)/d2*nphi*tan(pi/8.+phiav/4.)*tanphi 

par2=(2.*tanphi+2*nphil+nphi2)/(nphil*tanphi+nphi-l) 

par3=nphil*tanphi+nphi-l 

if  (c.le.0.0)  then 

cohesionless   soil 
P(i)=gamma*z/nphi*(dl*(dl/d2)**par3*exp((dl-d2)/d2*nphi* 

1   tanphi*tan(pi/8+phiav/4))-d2) 

else   if    (phiav.le.0.0)    then 

cohesive     soil 
P(i)=c*(dl*(3*alogl0(dl/d2)+(dl-d2)/d2*tan(pi/8)) 

1   -2*(dl-d2))+gamma*z*(dl-d2) 

else 

c-phi      soil 
P(i)=c*dl*(dl/d2)**par3*(l/(nphi*tanphi)*(exp(parl)-2*nphil 

1  *tanphi-D+par2)-c*(dl*par2-2*d2*nphi2)+gamma*z/nphi* 

2  (dl*(dl/d2)**par3*exp(parl)-d2) 


c 
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endif 
z"ce 
10    continue 

return 
end 


DATA 


10. 

500. 

125. 

16.5  40. 

10. 

30 

3. 

7.5 

2.e+09 

5.e+04  5 

.e+04 

0. 

0. 

3  1 

3 
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30. 

4.e+08      4.e+08 
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RESULTS 

SOIL  PARAMETERS 

Above  The  Critical  Surface 

Friction  Angle  =  10.0  degrees 
Cohesion  =  500.0 
Unit  Weight  =   125.0 

Below  The  Critical  Surface 

Coefficients  of  Subgrade  Reaction 

ksO  =  0.500e+O5 

ksl  =  0.500e+05 

krO  =  0.400e+09 

krl  =  0.400e+09 

kopt  =  3 


Variation  of  the  subgrade  reaction  coefficients 

nl  =  0. 
n2  =  0. 


PILE  DIMENSIONS 

Length  of  pile  above  the  critical  surface  =  16.5 
Length  below  the  critical  surface  =  40.0 

Distance  between  the  critical  surface  and  the 

bedrock  (or  hard  soil  layer)  =  10.0 

Length  of  pile  embedded  in  the  bedrock  ■  30.0 

Pile  Diameter  =  3.0 
Center-to-Center  Distance  =  7.5 
Stiffness  =  0.200e+10 

Number  of  points  that  the  pile  length  below  the 
critical  surface  is  discretized  to  =  30 


Boundary  Condition — Hinged  Head 

Force  Diagram 

force  on  the  pile  head  =   4845. 

force  acting  on  the  pile  at  the  critical 

surface  =  14556. 
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z 

DISPL(Y) 

SLOPE(S) 

MOMENT (M) 

SHEAR(V) 

■16.5 

0.439e-06 

0.191e-02 

-0.300e+01 

0.827e-K)5 

■13.8 

0.511e-02 

0.176e-02 

0.207e+06 

0.671e+05 

11.0 

0.944e-02 

0.136e-02 

0.365e+06 

0.471e+05 

-8.3 

0.124e-01 

0.783e-03 

0.462e+06 

0.227e+05 

-5.5 

0.137e-01 

0.122e-03 

0.486e+06 

-0.621e+04 

-2.8 

0.131e-01 

-0.514e-03 

0.424e+06 

-0.396e+05 

0. 

0.110e-01 

-0.999e-03 

0.264e+06 

-0.774e+O5 

1.3 

0.953e-02 

-0.114e-02 

0.162e+06 

-0.753e+05 

2.7 

0.794e-02 

-0.122e-02 

0.633e+05 

-0.736e+05 

4.0 

0.629e-02 

-0.123e-02 

-0.337e+05 

-0.722e+05 

5.3 

0.467e-02 

-0.117e-02 

-0.129e+06 

-0.711e+05 

6.7 

0.317e-02 

-0.105e-02 

-0.223e+06 

-0.703e+O5 

8.0 

0.186e-02 

-0.875e-03 

-0.317e+06 

-0.698e+O5 

9.3 

0.833e-03 

-0.633e-03 

-0.409e+06 

-0.695e+05 

10.7 

0.172e-03 

-0.329e-03 

-0.502e+06 

0.679e+05 

12.0 

-0.439e-04 

-0.855e-04 

-0.228e+06 

0.170e+06 

13.3 

-0.564e-04 

0.668e-05 

-0.483e+05 

0.898e+05 

14.7 

-0.261e-04 

0.190e-04 

0.112e+05 

0.238e+05 

16.0 

-0.566e-05 

0.102e-04 

0.152e+05 

-0.154e+04 

17.3 

0.120e-05 

0.276e-05 

0.713e+04 

-0.511e+04 

18.7 

0.171e-05 

-0.143e-06 

0.159e+04 

-0.279e+04 

20.0 

0.814e-06 

-0.572e-06 

-0.303e+03 

-0.769e+03 

21.3 

0.186e-06 

-0.317e-06 

-0.461e+03 

0.303e+O2 

22.7 

-0.319e-07 

-0.891e-07 

-0.223e+O3 

0.153e+03 

24.0 

-0.518e-07 

0.246e-08 

-0.521e+02 

0.865e+O2 

25.3 

-0.254e-07 

0.171e-07 

0.800e+01 

0.248e+02 

26.7 

-0.607e-08 

0.983e-08 

0.139e+02 

-0.400e+00 

28.0 

0.830e-09 

0.286e-08 

0.694e+01 

-0.459e+01 

29.3 

0.156e-08 

-0.148e-10 

0.170e+01 

-0.268e+01 

30.7 

0.790e-09 

-0.513e-09 

-0.204e+00 

-O.795e+O0 

32.0 

0.198e-09 

-0.304e-09 

-0.421e+00 

-0.430e-02 

33.3 

-0.207e-10 

-0.918e-10 

-0.216e+00 

0.137e+00 

34.7 

-0.471e-10 

-0.143e-ll 

-0.551e-01 

0.829e-01 

36.0 

-0.246e-10 

0.152e-10 

0.517e-02 

0.256e-01 

37.3 

-0.657e-ll 

0.914e-ll 

0.131e-01 

0.670e-03 

38.7 

-0.192e-12 

0.246e-ll 

0.695e-02 

-0.474e-02 

40.0 

0.   e+00 

0 .   e+00 

0.433e-03 

-0.489e-02 

z 

g 

j 

o 
q  I 


